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Abstract 


This is the first one of a series of papers about zeta regularization of the 
divergences appearing in the vacuum expectation value (VEV) of several local 
and global observables in quantum field theory. More precisely we consider a 
quantized, neutral scalar field on a domain in any spatial dimension, with ar¬ 
bitrary boundary conditions and, possibly, in presence of an external classical 
potential. We analyze, in particular, the VEV of the stress-energy tensor, the 
corresponding boundary forces and the total energy, thus taking into account 
both local and global aspects of the Casimir effect. In comparison with the 
wide existing literature on these subjects, we try to develop a more system¬ 
atic approach, allowing to treat specific configurations by mere application of 
a general machinery. The present Part I is mainly devoted to setting up this 
general framework; at the end of the paper, this is exemplified in a very sim¬ 
ple case. In Parts II, III and IV we will consider more engaging applications, 
indicated in the Introduction of the present work. 
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1 Introduction 


As well known, zeta regnlarization treats divergent qnantities in qnantnm field the¬ 
ory introdncing a complex parameter, with the role of a regnlator, and defining the 
renormalized observables in terms of analytic continnation with respect to the regn¬ 
lator. This constrnction was first proposed by Dowker and Critchley [30], Hawking 
[58] and Wald [101] to renormalize local observables, snch as the vacnum expectation 
valne (VEV) of the stress-energy tensor; the nltimate pnrpose was the semiclassical 
treatment of qnantnm effects in general relativity (e.g., nsing the stress-energy VEV 
as a sonrce in Einstein’s eqnations). After the previonsly cited pioneers, a nnmber 
of anthors championed the zeta approach to treat local observables; let ns mention, 
in particnlar, Cognola, Zerbini, Elizalde [25l[26], and Moretti [TTll60] 17 ^ 17611771 175] 
who worked in cnrved spacetimes, and Actor, Svaiter et ah EEi laiHni E7] who 
worked on spatial domains with bonndaries in flat spacetime. 

The zeta strategy can be as well applied to global observables, snch as the VEV of 
the total energy; in this version, it has in fact become more popnlar than its local 
connterpart. The literatnre on global zeta regnlarization is enormons; here we only 
cite the classical papers [88l[89] by Zimerman et ah, [13] by Blan, Visser, and Wipf 
and the monographies of Elizalde et al. [miMlIM], of Bor dag et al. [sms] and of 
Kirsten [62] (see also [211 mi SH EZl ESI E2]). 

Both in the local and in the global version, zeta methods provide a very natnral 
approach to stndy the effects of qnantnm vacnnm; in comparison with the orig¬ 
inal derivation of these effects by Casimir [20], and with other methods snch as 
point-splitting [T21 HSl [72] (in particnlar, see [55l [76] for a comparison between 
point-splitting and the zeta regularization approach) and the algebraic, microlocal 
approach (see [SI EH SSI EB Eg and citations therein), the zeta strategy is compet¬ 
itive and, perhaps, more elegant. 

The present series of papers (formed by this work and by the subsequent Parts 
II,III,IV [38l [39l SO]) considers both the local and the global zeta approach, devel¬ 
oping the viewpoint proposed in a special case by our previous work [36]. We are 
mainly interested in the theory of a quantized neutral scalar field on flat Minkowski 
spacetime, of arbitrary dimension d + 1. The held is conhned within a d-dimensional 
spatial domain H and arbitrary boundary conditions are prescribed for it on the 
boundary cAI; besides, we also admit the presence of a classical external potential, 
which could be meant to describe, in some sense, an effective theory of interacting 
helds. (We will only occasionally mention the possibility of replacing the hat domain 
n with a Riemannian manifold, or an open subset of it with prescribed boundary 
conditions; this amounts to replace Minkowski spacetime with a non-hat ultrastatic 
spacetime Ha-) 

Our attention is mainly focused on the VEV of the stress-energy tensor, of which 
we consider both the conformal and the non-conformal parts; however, we also 
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determine the total energy and boundary forces. 

Most of the works on local zeta regularization cited above consider Euclidean quan¬ 
tum fields; this makes a difference with respect to our formalism, based on canonical 
quantization in a genuinely Lorentzian framework. Moreover, each one of the pre¬ 
viously mentioned works [H El SI EHl ESI EZ] about local zeta regularization for 
flat domains deals with a specific spatial configuration (e.g. a strip between paral¬ 
lel planes, conhgurations with perpendicular planes, a rectangular wave guide or a 
rectangular box). In the present series of papers we are trying to develop a more 
systematic approach to the zeta strategy, to be applied nearly automatically in each 
specific case. 

In order to set up the general formalism we use with more generality the following 
fact, emerging from the previous literature: the analytic continuations involved 
in the zeta approach are closely related to certain integral kernels determined by 
the basic elliptic operator A which governs the spatial part of the field equation; 
among these kernels one could mention, in particular, the Dirichlet and heat kernels 
corresponding, respectively, to complex powers of A and to the exponential e~^. 
Our papers I-IV emphasize as far as possible the basic role of these and other kernels 
in view of zeta regularization. Our aim is to write down few very general rules to 
construct the required analytic continuation via integral kernels; these prescriptions 
can be applied in an almost mechanical way to treat specific configurations, as shown 
by several applications. 

The present Part I is mainly devoted to the general formulation, and in particular 
to the set of rules mentioned before; at the end of this paper we discuss, as a first 
application, the very simple case of a massless field on a segment (i.e., in spatial 
dimension d = 1) for several types of boundary conditions. In the subsequent 
Part II [5B] we will show how our general schemes work in a number of explicitly 
solvable cases (the field between parallel or perpendicular planes, or inside a wedge, 
with arbitrary boundary conditions and arbitrary choices of the spatial dimension 
and of the conformal parameter). In Parts III |39] and IV [iQ] we will consider 
two cases in which the implementation of our general rules requires a mixture of 
analytic and numerical calculations; more precisely. Part III will discuss a held in 
precence of a background harmonic potential and Part IV a held conhned within 
a rectangular box. In all these cases we will illustrate connections with the past 
literature including, when they exist, previous computations based on alternative 
approaches such as point-splitting. 

Let us describe in more detail the contents of the present Part I. In Section [2] 
and in the related Appendix we introduce the general framework for a neutral 
scalar held on a d-dimensional spatial domain hi, including a brief discussion on 
the stress-energy tensor with its conformal and non-conformal parts. This is an 
occasion to hx the attention on the basic elliptic operator A := —A -|- V(x) on 
n, where A is the Laplacian and V(x) the external potential. Always in Section 
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[21 we introduce zeta regularization in the formulation already used in our previous 
work [36]; the basic idea is to replace the quantized held 0 with its regularized 
version 0“ := where m G C is the regulator and k > 0 is a mass 

parameter. Formally, 0“ becomes 0 for m = 0; the zeta approach implements this 
idea in terms of analytic continuation. This means the following: for any local or 
global observable (say, the VEV of either the stress-energy tensor or of the total 
energy), after introducing a regularized version of the observable based on 0“, we 
dehne the renormalized value as the analytic continuation at k = 0. We distinguish 
between two versions of this prescription: the restricted zeta approach, in which 
the observable is analytic at m = 0, and the extended approach, where a singularity 
appears at m = 0 and is eliminated removing the negative powers of u in the Laurent 
expansion. 

In Section [31 (and in the related Appendices E] O and [E|) we present a number of 
integral kernels associated to the operator A := —A -|- ld(x). The Dirichlet kernel 
T)s(x,y) := ((5x|^“®5y) (s in a complex domain, x, y G ff) is closely related to the 
stress-energy VEV. More precisely, the regularized stress-energy VEV, built from 
0“, is determined by the Dirichlet kernel ^^(x, y) (and by its spatial derivatives) 
at y = X, s = (m ± l)/2; so, the renormalized version of this VEV is determined 
by the analytic continuation of Ds near s = ±1/2. As shown in the cited section, 
the continuation of the Dirichlet kernel can be determined algorithmically relating 
it to the heat kernel iF(t;x, y) := {5:x\e~^5y) or to other kernels (among them the 
so-called cylinder kernel T(t;x, y) := ((5x|e“‘'^(5y), considered by Fulling [331 [37]). 
This precedure ultimately gives the set of mechanical rules for zeta regularization 
mentioned previously in this Introduction. 

In Section [Hand in Appendix [F] we relate to the previous framework the total energy, 
the pressure and the total forces on the boundary. We also take the occasion to prove 
(at the regularized level) the equivalence between two alternative dehnitions of the 
pressure, often assumed without proof in the literature: pressure as the action of the 
stress-energy VEV on the normal to the boundary, and pressure as the functional 
derivative of the bulk energy with respect to deformations of the spatial domain D. 
In Section [5] and in Appendix [G] we consider some variations of the general frame¬ 
work of Sections [2] and [3] concerning the following situations: i) the case where 0 
is either an isolated or non-isolated point of the spectrum of the fundamental op¬ 
erator A := —A ± E(x); ii) the case where the flat spatial domain D is described 
via curvilinear coordinates or, more generally, the case where D is a (possibly non¬ 
flat) Riemannian manifold equipped with arbitrary coordinates (or an open subset 
of it, with prescribed boundary conditions). This suggests the possibility to apply 
our formalism to the case of a non-flat ultrastatic spacetime where the line element 
reads ds^ = —dt^ ± with the Riemannian line element of D. 

The hnal Section [6] presents a hrst application of our formalism; this concerns the 
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case in which d = 1 and is the interval (0,a), with suitable boundary condi¬ 
tions. This conhguration is very simple: we take it in consideration just to show in 
action our mechanical procedures for analytic continuation. Comparison with the 
existing literature on this simple case is performed. As already mentioned, more 
sophisticated applications of our approach will appear in Parts II-IV. 

2 Zeta regularization for a scalar field 

2.1 General setting. In this section we summarize the zeta regularization method 
for the propagator and the vacuum expectation value (VEV) of the stress-energy 
tensor of a quantized scalar held, in the formulation presented in [SB] (see also [BT]). 
Here the scheme of is slighlty generalized, admitting the presence of a classical 
background potential V and arbitrary spacetime dimensions. 

Throughout the paper we use natural units, so that 

c=l, n=l. (2.1) 

We work in (d -j- l)-dimensional Minkowski spacetime; this is identihed with 
using a set of inertial coordinates 

X = = (x°,x) = (t,x) , (2.2) 

in which the Minkowski metric has coefficients {rjnu) = diag(—1,1,..., 1). Let us 
hx a spatial domain C where we consider a quantized neutral, scalar held 0 
in presence of a classical background static potential V; so, we have H R, 

X I—)■ 1/ (x) and 

0 : R X £^(5") ; 0 = {-dtt A - 1/(x))0(x, t) , (2.3) 

(analogous settings are considered, e.g., in [HI [151HTJ 06] ) (15 • Here we are referring 
to the space /I(5^) of linear operators on the Fock space and to the subset Csa{d) 
of selfadjoint operators; A := d-dimensional Laplacian. We assume 

V to be smooth, and prescribe appropriate boundary conditions (e.g., the Dirichlet 
conditions 0(t, x) = 0 for x G dfl). From here to the end of the paper we put 

A:=-A + V, (2.4) 

intending that the boundary conditions are accounted for in the above dehnition; 
we assume this framework to grant that M is a selfadjoiont operator in the Hilbert 

^ Of course the notation 0 : R x P —>• TsalS^), (x, t) !->■ 0(x, t) is used here in connection with a 
generalized operator valued function; in fact, as well known, 0 is an operator valued distribution. 
In the papers of this series we use the classical functional notations even for generalized functions. 
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space with the inner product 


(fig) ■= [ dxf{x)g{x) (2.5) 

Jn 

(dx denoting the standard Lebesgue measure on f2) ([^). Moreover, we assume A 
to be strictly positive, by which we mean that the spectrum (t{A) is contained in 
[£^,+cxo) for some £ > 0. Let us mention that A may have continuous spectrum, 
a fact typically occurring when is unbounded; therefore, when we speak of the 
eigenvectors of A we always intend them in a generalized sense, including improper 
eigenfunctions. 

Let us consider a complete ortho normal set {Fk)k^ic of (generalized) eigenfunctions 
of M (0), indexed by an unspecihed set of labels /C, and let us write the corresponding 
eigenvalues in the form {ojDkeic ^ ^ for all fc G /C). Thus 

Ffc : —)■ C; AFk = io^Fk ; 

{Fk\Fh) = 5{k, h) for all fc, h G /C . ^ ‘ ^ 

Any eigenfunction label k E K, can include different parameters, both discrete and 
continuous. Besides, we generically write dk to indicate summation over all labels 
(i.e., literal summation over discrete parameters and integration over continuous 
parameters, with respect to a suitable measure); 6{h, k) = 6{k, h) is the Dirac delta 
function for the label space /C (this reduces to the Kronecker symbol in the case of 
discrete parameters). Note that the condition cjfc ^ e > 0 excludes the presence of 
infrared divergences from all the sums over k appearing in the sequel. 

The functions 


/fc : R X D C , X = (f, x) H- /fc(a;) := e *‘^''*Ffc(x) (2.7) 


fulhll {—dtt — A)fk = 0 ; they allow us to infer for the quantized held a normal modes 
expansion of the form 


0(x) 



akfk{x) + a\fk{x) 


( 2 . 8 ) 


^In passing, we recall that the Hilbert space "S can be realized as the direct sum of all sym¬ 
metrized tensor powers of L^(H). 

^For a fully rigorous discussion of generalized eigenfunctions, see Chapter IV of [S2]. In the 
sequel, following the usual terminology, when speaking of functions (or distributions) on ti we 
use the adjectives “proper” or “improper”, to distinguish between objects which actually belong 
to FitX) or not. In this spirit we speak of proper or improper eigenfunctions, and use the same 
terminology for the corresponding eigenvalues. In the sequel, the adjective “generalized” in relation 
to eigenfunctions is sometimes omitted. 
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(with ^ indicating the adjoint operator, and the complex conjugate). In the above 
we are considering the destruction and creation operators a/c, a|, 6 which fulhll 

the canonical commutation relations 

[ak,ah\ = 0 , [afc,a[] = 6{h,k) , a^lO) = 0 , (2.9) 

where |0) G 5^ is the vacuum state (of unit norm). A relevant character for the sequel 
of our analysis will be the propagator, i.e., the vacuum expectation value (VEV) 

(O|0(a:)0(2/)|O) {x,yeRxn). (2.10) 

Let us pass to the stress-energy tensor operator; this depends on a parameter G R, 
and its components T^i, : R x hi —£ 50 ( 5 ^), for /i, G {0,1,..., d}, are given by 

7 ;, := (1 - 20 o dA - Q - 2^ +V^)-2i^o (2.11) 

where we use the symmetrized operator product Ao B := (1/2) (AR -|- BA) and all 
the bilinear terms in the field are evaluated on the diagonal (e.g., indicates 

the map x 1 — )■ d^(j){x) o dy4>{x)). 

Eq. (12.lip provides a natural quantization of what is often called the “improved” 
stress-energy tensor; this is a well-known modification of the canonical stress-energy 
tensor with an additive term proportional to the parameter 0 fhat does not alter 
its divergence. For certain boundary conditions, e.g. of Dirichlet type, this addition 
does not even alter the corresponding momentum vector. We refer to Appendix 
El for further details on this topic. Here we only recall that the improved stress- 
energy tensor was hrst proposed by Callan, Coleman and Jackiw [19] in order to 
deal with some pathologies appearing in perturbation theory; later on, this tensor 
was reinterpreted in terms of the Minkowskian limit for a scalar held coupled to 
gravity via the curvature scalar Ha El EQl Ea El]. 

Needless to say, the vacuum expectation value of the stress-energy tensor (12.lip is 
the main character in the theory of the Casimir effect. It is evident from Eq. (12.lip 
that (0|T^i,(x)|0) can be expressed formally in terms of the propagator (I2.10p (and of 
its derivatives) evaluated on the diagonal y = x. On the other hand, the propagator 
is known to be plagued with ultraviolet divergences along the diagonal, so that 
(0|7/iy(a;)|0) is a merely formal expression for a divergent quantity; our purpose is 
to redehne the propagator and the stress-energy VEV via a suitable regularization, 
ultimately yielding hnite values for these quantities. 

2.2 Zeta regularization. Let k > 0 denote a parameter, to which we attribute 
the dimension of a mass (or of an inverse length, since h = 1). k will be called 
the mass scale; it is introduced for dimensional reasons and plays the role of a 
normalization scale. See Ha EH Ea EDI [75] for further comments regarding this 
parameter and its presence or absence in the renormalized observables of the field. 










We will check that the hnal, renormalized results depend on k only when singularities 
appear in the analytic continuations involved in the following construction. 

The zeta strategy, in the version proposed in [36] to give meaning to the VEV of 
relies on the powers 

( 2 . 12 ) 

where A is the operator fl2.4l) and m G C ; these are employed to dehne the smeared, 
or zeta-regulahzed, Geld operator 

0 “ ;= 0 , (2.13) 


depending on the complex parameter u and coinciding with the usual held operator 
(j) for M = 0. If (Tfc)fcgyc is a complete orthonormal set of eigenfunctions of A with 
eigenvalues {uDk&fc, we have {K~‘^A)~'^^^Fk = for any k & K.] so, 

the functions /^(f,x) = fulhll {K~'^A)~'^^‘^fk = and the 

expansion fl2.8p for 0 becomes, after application of {k~^A)~^^‘^, 


0“(a;) = 




dk 


Ik 


^/2 


u 


1/2+U/2 


dk , 


+ 


fk{x) 


(2.14) 


Note that, in the limit ook —)■ +cxo, the term 1 / 0 ;^^''"“^^ in the above integral vanishes 
rapidly if 3?^ is large; this is a manifestation of the regularizing effect of the operator 
for large a fact we will describe much more precisely in the sequel. 
Using 0“, we can dehne a regularized propagator 

(O|0“(a;)0“(|/)|O) {x,yeRxn) (2.15) 


and a zeta regularized stress-energy tensor 

yi:= (i-209y'‘oaJ'‘-f^-2ij,,^,(ay'‘a,r+v(rfA2(rod,J'‘, (2.16) 

where, as in Eq. fl2.1ip . all the bilinear terms in the held are evaluated on the 
diagonal. 

We are interested in the VEV of this regularized stress-energy tensor, which formally 
gives (0|T^i,(x)|0) in the limit m —)■ 0. Of course, we can relate the VEV of T“^,(a;) 
to the regularized propagator (12.151) in the following way: 

(o|f;t(i)|o) = 

^2 ^ (dx9y‘'-h ^2 cIj^A-|-U(x)^ dya-yv) 

• (O|0“(x)0“(2/)|O) . (2.17) 
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We will return later on this equation and on its use for the actual computation of 
the above VEV. 

Typically, the regularized propagator and the VEV of are analytic functions 

of M, for sufficiently large; the same can be said of many related observables 
(including global object, such as the total energy, which is related to the space 
integral of the (0,0) component of the stress-energy tensor). Let us consider any 
one of these (local or global) observables, and denote with J^{u) its zeta-regularized 
version, based on Eq. fl2.13p (see, e.g., Eq. fl2.15p or Eq. fl2.16p i: we assume that 
the function u ^ F{u) is well dehned and analytic for u in a suitable domain Uq 
of the complex plane. The zeta approach to renormalization can be formulated in 
either a “restricted” or an “extended” version, both described hereafter. 

i) Zeta approach, restricted version. Assume that the function Uq ^ C,u J^{u) 
can be analytically continued to a larger open subset W of C such that 0 G W ; let us 
use the notation u ^(u) even for this extension (j^). In this case, making reference 
to the analytic continuation, we define the renormalized value of the observable 
under consideration as 

J^ren := -T(0) . (2.18) 

ii) Zeta approach, extended version. Assume that there is an open subset W of C, 

larger than Wq, such that 0 G W and the function u G Uq iF(u) has an analytic 
continuation to W\ {0}, still indicated with J^. In this case, since there is an isolated 
singularity at u = 0, in a neighborhood of this point we have the Laurent expansion 
J^(u) = ■ Let us consider the regular part 

+ 0O 

(/2F^)(u) ; (2.19) 

k=0 

we define the renormalized value of the given observable as 

J^ren := (FPF)(0) (2.20) 

(i.e., Tren = -Tq)- In most applications F is meromorphic close to u = 0, which 
means that it has a pole at this point; in this case the previous Laurent expansion 
has the form J^{u) = i where N G {1, 2, 3,...} is the order of the pole. 

Let us stress that the prescription fl2.20p is a quite straightforward generalization of 
the approach considered in [T31 [33], where T was assumed to posses a simple pole 
in u = 0 (i.e., it was assumed that N = 1). 

®In the style of our previous work [36] we should write iF : Uq ^ C for the initially given 
function and AC jF : U ^ C for its analytic continuation; here, we choose to simplify the notation, 
writing if for both functions. Note that the analogue of Eq. (12.181) in the style of [3B] would be 

Tren := (dC J-)(0). 
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Of course, the restricted zeta approach of item (i) is equivalent to a special case of 
the extended approach, in which J^{u) has a removable singularity at u = 0 and the 
Laurent expansion at this point is the usual power series expansion. 


A large part of our subsequent work will be devoted to the application of the pre¬ 
vious scheme to the VEV of the stress-energy tensor. In general, we define the 
renormalized version of the latter as 


{Q\T,,{xmren-.= RP (0|T“(x)|0) ; 


u=0 


( 2 . 21 ) 


when no singularity appears at u = 0, according to the restricted approach (i) the 
above definition reduces to 


(0|j;„(x)|0)„„ := (0|f“ (i)|0) 


( 2 . 22 ) 


In [32], we only considered the prescription fl2.22p in the special case of a Dirichlet 
field (with 1/ = 0) between two parallel planes, i.e., in the configuration correspond¬ 
ing to the standard theory of the Casimir effect. In that case the approach fl2.22p 
was implemented via a direct computation of the analytic continuation appearing 
therein; as already stressed, here we aim to much more generality. 


2.3 A remark. In the sequel, while performing zeta regularization and the con¬ 
sequent renormalization, it is sometimes natural to consider, in place of u, some 
complex parameter s related to n by a simple transformation. In view of such sit¬ 
uations, it is convenient to generalize some notations of the previous subsection in 
the following way: 

i) Consider an analytic function iSq —)■ C, s i—)■ J^(s), where So is an open subset of 
C; if this admits an analytic continuation to a larger open subset S, the latter will 
be still denoted with s i—)■ J^(s). 

ii) Suppose the analytic function iSq —)■ C, s i—)■ J^(s) has an analytic extension to 

S \ {so}, where S is an open subset of C and Sq E S . Then, the Laurent expansion 
R{s) = '^k=-oo^kis — So)’^ will be used to define the regular part (near sq) of this 
analytic continuation as {RPR){s) := .Tfc(s — sq)^ ; of course, this implies 

{RPR)iso)=Ro- 

2.4 Staticity features of the VEV of Let us return to the regularized 

stress-energy tensor; for x = (f, x) G R x we claim that 

{•)|tk(x)|0) is independent of ( , 

(0rax)|0) = {0|i;“„(i)|0)=0 fort 

These statements are not surprising, due to the staticity of the general framework 
considered in the present paper; a formal proof will be given in subsection 13.81 (see 
Eq.s fl3.31113.33P and the considerations which follow them). Of course the features 
of Eq. fl2.23p are preserved by analytic continuation, so that an analogue of this 
equation holds for the renormalized VEV (0|T^,o(x)|0)ren as well. 
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2.5 Conformal and non-conformal parts of the stress-energy tensor. In 

the literature (see, e.g., m [m M) it is customary to write the stress-energy 
tensor (here to be intended as one of the operators or either one of the 

VEVs {0|T“y|0), (0|T^i,|0)ren) as the sum of a conformal and a non-conformal part. 
In order to dehne these quantities, let us consider for ^ the critical value 


■= 


d — 1 


(2.24) 


it is known that, when coupling of the scalar held to gravity is taken into account, 
the theory is invariant (for C = 0) under conformal transformations of the spacetime 
line element if ^ has the above critical value (see, e.g., |102] . p.447). In the sequel 
we adopt the notations 


0 = conformal , ■ = non-conformal 


(2.25) 


and we dehne, for example, the conformal and non-conformal parts of the renormal¬ 
ized stress-energy VEV (0|7)ji.|0)ren in the following way: 


(0|fi^)|0),en:=(0|f^.|0),en 

{0|f™|0)„„ := ^ ({0|f,„|0)„„ - (0|f^W|0)„„) . 

Of course, this implies 

(0|f^.|0),en = (0|fW)|0),en + (^-^d) (0|fW|0)_ . 


(2.26) 

(2.27) 

(2.28) 


In the applications to be considered in Section [6] and in the subsequent Parts II,III 
and IV, when presenting our hnal results for the renormalized stress-energy VEV, 
we will either write them in the form (I2.28p or give separately the conformal and 
non-conformal parts (I2.26p (I2.27p . 

2.6 Total energy and pressure on the boundary. The total energy is, by 
dehnition, the integral of (0|Too(a:)|0) over the spatial domain O. We defer the 
discussion of this topic to Section 01 therein we will describe the representation of 
the total energy as the sum of a bulk term and a boundary term, in the framework 
of zeta regularization. 

In the same section we will use zeta regularization to treat the pressure on the 
boundary dfl of the spatial domain; this quantity can be dehned in terms of the 
VEV of the spatial components Tjj . There is an alternative characterization of 
the pressure in terms of the variation of the bulk energy (see Eq. fl4.6p for the 
dehnition) with respect to deformations of the spatial domain ; the equivalence 
of this dehnition with the previous one has often been assumed uncritically in the 
literature, so we think it can be useful to produce a formal proof (see Section 0]). 
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3 Expressions of the zeta regularized stress-energy 
VEV in terms of integral kernels 

In this section hi always denotes a spatial domain in and we often consider the 
Hilbert space LP‘{VL) of the square integrable complex-valued functions on hi. 

3.2 Basics on integral kernels. Let us consider a linear operator B acting on 
The integral kernel of B is the (generalized) function 

B{,):nxn^C, (x,y)^H(x,y):=(4|S^y) (3.1) 

where (5x and 6y are the Dirac delta functions centered at x and y, respectively, 
here viewed as improper vectors of the Hilbert space (Q). Equivalently, the 

integral kernel of the operator B can be defined as the unique (generalized) function 
B{ , ) : D X D —)■ C such that 

{Bij){x) = [ dy H(x,y) V'(y) , (3.2) 

Jn 

for all sufficiently regular ip : Q ^ C. \i B possesses a complete orthonormal 
set of (generalized) eigenfunctions {Fk)keK. with corresponding eigenvalues /3k & C 

®Let us repeat what we have declared in the footnote of pageO in these papers we use functional 
notations even for objects which are not ordinary functions. Making reference to the theory of 
Schwartz distributions, we can understand the statement “B{ , ) : 17 x H —>■ C” as a way to indicate 
that we are considering a (complex) distribution on 17 x 17; let us describe how to intend Eq. (13.11) 
from this viewpoint. We refer to the space of test functions 72(17), formed by the C°°, compactly 
supported functions : 17 —>■ C and equipped with its standard inductive limit topology; 17(17 x 17) 
has a similar meaning. We wish to define rigorously (x, y) >->• {Sx\B Sy) as a distribution on 17 x 17, 
i.e., as a continuous linear form on 77(17 x 17). To this purpose we consider y, b S 77(17) and note 
that the identities ‘P = Jq dx(p(x)Sx and ip = Jq dx'i/(x)Sx give formally 

[ dxdy if{x){6x\B6y)tp{y) = (iplBip) . 

On the other hand, if the domain of the operator B contains D{^1) and the sesquilinear map 
77(17) X 77(17) —>• C, {ip, (p) >->• {(p\Bip) is continuous, using the nuclear theorem of Schwartz (see 
|92) . Thin. II) one proves rigorously the existence of a unique distribution on 17 x 17, denoted with 
(x,y) !->• (5x|-B(5y), such that the above relation holds for all Pjtp £ 77(17), intending the integral 
in the left hand side as the action of this distribution on the test function (x, y) i—>■ Tp{x)ip(y). 
One can use similar considerations to give a distributional meaning to Eq. (13.21) (see again [^1. 
Throughout the papers of this series, whenever we speak of the integral kernel of an operator 
B we implicitly assume B to possess the regularity features mentioned before, so that {6x\B6y) 
makes sense at least as a distribution. In some cases described explicitly in the sequel, stronger 
assumptions on B ensure that {Sx\BSy) can be understood as an ordinary function, continuous with 
its derivatives up to a certain order; this is the situation outlined in Appendix [B] often mentioned 
in the sequel in relation to several kernels of interest for us. 
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{BFk = /3kFk), then 

;B(x,y)= [ dkl3kFk{x.)Fk{y) (3.3) 

Jk 

(since the function in the right-hand side fulhlls equation fl3.2p for all 'ip). The 
precise sense in which the eigenfunction expansion fl3.3p converges depends on the 
specihc features of the operator B] we generally assume distributional convergence 
(P), leaving to future subsections the consideration of special cases where convergence 
can be intended in a stronger sense. 

Incidentally, let us mention the relation existing between the kernel B( , ) and the 
trace of B; the latter, if it exists, is the number TiB := Jj^dk {Fk\BFk) G C, where 
{Fk)k&K is any complete orthonormal set of LP‘{VL). The right-hand side does not 
depend on the choice of {Fk)keK.', in particular, if B has purely discrete spectrum, 
(Fk) k^fc is a complete orthormal set of proper eigenfunctions labelled by a countable 
set fC and {(3k)keK are the corresponding eigenvalues, we have TtB = Yhk&Kl^k-i h 
this series converges. Returning to dehnition fl3.ip of the kernel B{ , ), we see that 

Tri3= f dxB{x,x) (3.4) 

Jn 

(since (5x)xeo is a generalized complete orthonormal set) (P). 

Let us move on and note that the boundary conditions possibly involved in the 
dehnition of B have implications for the kernel B{ , ); for example, if boundary 
conditions of the Dirichlet type are involved, the eigenfunctions (T). )fc6/c in Eq. ([23]) 
vanish on dfl, thus yielding i3(x, y) = 0 for x G dQ or y G dfl. 

Let us also mention that from Eq. fl3.ip one infers 

Si(x, y) = B{y, x) , S(x, y) = S(x, y) (3.5) 

where B^ is the adjoint operator of B with respect to the inner product of L^(ff) 
while B is the complex conjugate operator, such that Bpj = B ip for all ip. These 
facts imply 

B{y,x) = B{x,y) if B^ = B . (3.6) 

3.3 The operator .4.. Most of the integral kernels considered in the sequel will 
be related to a selfadjoint operator A acting in L^(f2); for example, they will be the 
kernels associated to some function of the operator A (say, a power M^). To treat 
these kernels, precise assumptions on A will be specihed whenever necessary; in any 

^Meaning that J^^^dxdy B{x,y)A^,y) = J^dkpk J^^^dxdy Fk{x)Fk{y)if{x,y) for all test 
functions tp on x Q. 

®In the sequel, the adjective “generalized” in relation to complete orthonormal sets is sometimes 
omitted. 
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case, we will typically consider three situations of descreasing generality, described 
by the forthcoming Eq.s fl3.7pfl3.8pfl3.9l) . 

In the hrst situation, we will simply assume that 


^ is a strictly positive, selfadjoint operator in L^(f2) (3.7) 

(let us recall that strict positivity means cr^A) C +cxo) for some £ > 0). 

In the second situation, the assumptions are the following: 

A = —A + 1/ in with V a real C°° potential on and , . 

boundary conditions such that A be selfadjoint and stricly positive. ^ ' 

In the third situation, the assumptions are: 

A = —A + V on a. bounded domain with dfl of class and Dirichlet , . 

boundary conditions, V a real C°° potential on and V (x) ^ 0 for all x G ^ ' 


Here and in the sequel U dfl is the closure of The assumptions fl3.9p 

grant selfadjointness and strict positivity of A; moreover, they imply that A has 
a purely discrete spectrum and one can build a complete orthonormal system of 
proper eigenfunctions with eigenvalues ojI labelled by /C = {1, 2, 3,....} in such a 
way that 0 < Wi ^ a ;2 ^ tua ^ ... (with the possibility that some of these inequalities 
are equalities, to deal with the case of degenerate eigenvalues). It is well-known that 
the eigenvalues, when ordered in this manner, fulhll the Weyl asymptotic relation 

(Ufc ~ C for k —)■ -l-cx) , (3.10) 


where C := 2v^r((i/2-|-Vol(fI)“^/'^ (see [H], Thm.5, page 189 and [3^, § 
8.2, pages 99-101 for elementary derivations). Appendix iBl contains a number of 
technical results concerning cases fl3.7p fl3.8p fl3.9p . that will be mentioned whenever 
necessary in the sequel of this section. 

These results often refer to the spaces C-^(fl) and C^{fl x f2) or (in the case fl3.9p j 
C^{fl) and C^{fl x fl), for j G N or j = cxo; these are formed by the complex func¬ 
tions on fl, X and so on, which are continuous along with their partial derivatives 
of all orders ^ j. In particular, in the rest of the section we will present situations 
connected to cases fl3.7p fl3.8p fl3.9ll . in which certain integral kernels (x, y) i-G- i3(x, y) 
related to A are of class C^{fl x f2) or C^{fl x fl), for suitable j. 

The condition of strict positivity for A in Eq.s fl3.7p fl3.8p will be occasionally re¬ 
laxed, assuming only that the eigenvalues of A are non-negative (and declaring this 
explicitly, to avoid misunderstandings); see, e.g., subsection 13.101 

3.4 The Green function. Let .4, be a strictly positive selfadjoint operator in 
L^(f2); then, we can introduce the inverse operator A~^ and the corresponding 
kernel 

G(x,y) := ^■^(x,y) , (3.11) 
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which is called the Green function of ^. In terms of this kernel, the identity AA ^ = 
1 can be re-expressed as 

^xG'(x,y) = 5(x-y) , (3.12) 

where indicates the operator A acting on G(x, y) as a function of x. Using 
a complete orthonormal system {Fk)k&K of eigenfunctions of A with corresponding 
eigenvalues {uDkeic, we can express the Green function as 

r _ 

G'(x,y) = — Fk{^)Fk{y) . (3.13) 

Jk 

The Green function is among the most familiar integral kernels, especially when 
A = —A+V with suitable boundary conditions (of the Dirichlet, Neumann or 
Robin type); in this case, uniqueness results are available for the Poisson equation 
(perturbed with an external potential), allowing to characterize the Green function 
G(x, y) as the unique solution of Eq. fl3.12p fulhlling the prescribed boundary con¬ 
ditions for X G dfl or y G dfl. The literature on this topic is enormous, and here we 
only mention some well-known monographies: Berezanskii [9], Krylov [M], Sauvi- 
gny [90] and Shimakura [96] give abstract and rigorous analyses, while Duffy [3T] . 
Kythe [65], Sommerfeld [99] and Stakgold and Holst [100] present more practical 
and explicit discussions. 

3.5 A digression on complex powers. Throughout the present paper (and in 
Parts II-IV), the following conventions are employed: 

i) In : (0, -|-cxd) — )■ R is the elementary logarithm; 

ii) for any a G C, we systematically refer to the standard dehnition 

:= for all x G (0, +oo) ; (3.14) 

hi) for a G C and z in a convenient subset of the complex plane, we dehne 

^aln\z\+iaa,rgz /g 


where arg : —)■ R is some determination of the argument; this determination 

depends on the domain and must be specihed in each case of interest. In most 
applications considered hereafter we set 

CV= C \ [0,+<x.) i 

arg := the unique determination of the argument with values in (0, 27r) . 


3.6 The Dirichlet kernel. Let again A be a strictly positive selfadjoint operator 
in L^(f2). The power can be dehned through the standard functional calculus 
for each s G C; the corresponding integral kernel 

D,(x,y) := A“"(x,y) (3.17) 
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is called the s-th Dirichlet kernel. In passing, let us note that D_i{'x,y) coincides 
with the Green function G(x, y) considered in subsection 13.41 

If {Fk)keK is a complete orthonormal set of eigenfunctions of A with corresponding 
eigenvalues {coDkeic we have A~’^Fk = so that (by Eq. fl3.3p l 

r dk _ 

^s(x,y)= / — Ffc(x)Ffc(y) . (3.18) 

JK. 

The denomination of “Dirichlet kernel” employed for Dg is suggested by the simi¬ 
larity between the above expansion and the Dirichlet series, considered in [561 1731 

For our purposes, it is important to give sufficient conditions under which Ds(x, y) 
is a regular function of (x, y) (even on the diagonal y = x, of special interest in the 
sequel); we are also interested in cases in which the expansion fl3.18p converges in a 
stronger sense, in comparison with the distributional sense that we are generically 
ascribing to kernel expansions. Let us present some results of this kind, which are 
proved in Appendix |Bl 
Let A = —A -|- D be as in Eq. fl3.8p : then 

d 1 

Dg E C^iVL X D) for s e C, j G N such that > - + ^ . (3.19) 

With the stronger assumptions fl3.9p . we have 

Dg E C^iytx D) for s G C, j G N such that > — + ^ . (3.20) 

Again with the assumptions fl3.9|) . A possesses a complete orthonormal set of proper 
eigenfunctions {Fk)k=i,2,3,... (see the comments after the cited equation), and the 
expansion fl3.18p reads 


+ 00 

Ds{^,y) = 


1 , 


, ,2s 

k=l 


Ffc(x)Ffc(y) 


(3.21) 


In this case one has Fk G (^““(D) for each fc; moreover, for > d+j/2 the expansion 
fl3.2ip is absolutely and uniformly convergent on D x D, with all its derivatives up 
to order j\ see Appendix [B] (especially Eq. (1B.26P ) for more details. 


To conclude the present subsection let us mention that the general statement of Eq. 
(13.4p . here applied with B = A~^, yields 


dx D Ax, x) = Tt A 


(3.22) 
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provided that the above trace exists. With the assumptions fl3.9|) for A = —A + V, 
using the Weyl estimates fl3.10l) for cui ^ a ;2 ^ ••• we hnd 


+ 00 , , 

Tr^“^ = is hnite for 3fJs > - . (3.23) 


3.7 Some remarks concerning the Dirichlet kernel and its derivatives. 

Let us consider again a strictly positive selfadjoint operator A in moreover, 

assume this operator to be real, in the sense that .4, = .4, (i.e., Aip = A'lp for all ip). 
If {Fk) k£ic is a complete orthonormal set of eigenfunctions of A with related eigen¬ 
values {ul)k£K, then the conjugate system {Fk)k£K is as well a complete orthonormal 
set of eigenfunctions of A with the same eigenvalues; so, besides Eq. fl3.18p we have 
an alternative representation for the Dirichlet kernel Ds(x, y), based on this con¬ 
jugate system. From here, we easily infer that, for s G C with complex conjugate 

s, _ 

D,(x,y) = D,(y,x) and ^.(x, y) = ^^(x, y) . (3.24) 

To go on we claim that, for any pair of multi-indexes a, 13, 

a“a^D,(x,y)|y^^ = dfd“D,(x,y)|y^^ . (3.25) 

Indeed, due to the hrst identity in Eq. fl3.24l) . we have d“d^Ds(x, y) = d°d^Ds{y, x); 
when evaluating the right-hand side of this equality on the diagonal y = x, the 
variables can be relabeled to yield ^“^^^^(y,x)= dyd^Ds(x,y)|^_^, thus 
proving Eq. fl3.25|) . 

All the above results can be applied to the (real) operator A := —A -|- E(x); the 
symmetry properties outlined here for the corresponding Dirichlet kernel will be 
relevant in connection with the results of the next section. 

3.8 The regularized propagator and stress-energy VEV: connections with 
the Dirichlet kernel. Let us refer to the framework of the previous section, where 
the operator A = —A -|- E in L^(f2) has been considered in connection with a 
quantized scalar held. In the sequel x = {x^,x.),y = {y^,y) G R x D; if we use the 
expansion fl2.14p for the regularized held (jA in terms of creation and destruction 
operators we obtain for the regularized propagator the expression 


dkdh 


= K 


kxK 2{ukUJhy^^ 


(O|0“(a;)0“(|/)O) = 

akfkix)+alfkix) ahfhiy)+alfhiy) | 0 ) • 


(3.26) 
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This relation, along with the identities (0|afca/i|0) = (0|a|,a|j|0) = (0|a|,a/i|0) = 0 and 
{0|afca|^|0) = 5{k,h), gives 


(O|0“(x)0“(|/)|O) = K 
f dk 


dk 


2 cjfc 


l+ti 


fk{x)fk{y) = 


= K 


2 cjfc 


l-\-u 


Ffc(x)Ffc(y) e 


-iujk(x°-y°) 


(3.27) 


From here we can easily obtain the derivatives of the propagator; for example, for 
j G {1, ...d}, we have 


d^Oyj{0\^^{x)^^{y)\0) 



dk 

2ujk^ 


Fk{x.){dyiFk){y) 




(3.28) 


In particular, if we apply Eq.s (I3.27p (I3.28p with y = x and compare with the 
eigenfunction expansion (I3.18p of the Dirichlet kernel, we get 


(O|0“(a;)0“(2/)|O) = ^ T>2^(x,y) 

y=x L 2 


y=x 


a^Oyi(O|(/)“(a;)0“(|/)|O) = -iiF dyjDy.{^,y)\ 


y=x 


y=x 


(3.29) 

(3.30) 


One can express similarly all the derivatives with y = x appearing in Eq. fl2.17p 
for the regularized stress-energy VEV. In this way (and using as well the identity 
(I3.25P ) we obtain the following results, where qj, ^ are spatial indexes ranging in 


{0|r“(«,x)|0)=K' 


^ +d-Dajl (x, y) + (j + r (x))Ds+i (x, y) 


(0|T5.(i,x)|0) = {0|7™(«,x)|0) = 0, 


. (3.31) 


y=x 


(3.32) 


^To prove Eq. (13.321) . note that 


(0|fo“(t,x)|0) =-(0|fo“(t,x)|0) =-^ (dyjE>i(x,y) - d^j£)-(x,y)) 


y=x 


and that the last expression vanishes due to identity (13.251) . Besides, note that Eq. (I3.33P is 
equivalent to the more symmetric expression 


(0|T“(hx)|0) = (0|T“(hx)|0)=Ki 
'1 ^ 


(4 - (d"^\^-hF(x))£)^(x,y)) 

+ + - I + dyiyj^D^i{yi,y) 


19 






























(3.33) 


= 


(0|7;“(t,x)|0) = (0|i;“(t,x)|0) = 

^)^ii(.C)ji^(x,y) - {d'''dyi + V{x))D + 


~i~ ((2 ^^ (^)y) 


y=x 


The above equations indicate, amongst else, that (0|T“j,(t, x)|0) does not depend on 
the time variable t; this comes as no surprise at all, since our general framework is 
itself static (indeed, the spatial domain hi and the potential V are time indepen¬ 
dent). These features of the regularized stress-energy VEV had been anticipated in 
subsection 12.41 due to them, in the rest of the paper we will use the notation (j^ 


(0fc(x)|0)^(0fc(t,x)|0) . (3.34) 


Up to now we have been working rather loosely, but it is not difficult to indicate 
the precise conditions for the validity of our manipulations. Indeed, Eq.s (I3.31tl3.33|l 
for {0|T“j^(x)|0) involve the Dirichlet kernels and D^+i, with its second order 

derivatives, along the diagonal y = x. So, for these equations to be meaningful it 
is sufficient that E C^(Q x Q) and Du+i G C^(Q x Q); with the assumptions 

(13. 8 p on ^ = —A -1- V, recalling Eq. (I3.19p we see that both conditions are fulhlled 
if 


> d -|- 1 . 


(3.35) 


For A, u as in Eq.s (13. 8 p (I3.35p . the map x (0|T“^,(x)|0) is in C°(r2) for all /i, v. 
With the stronger assumptions (13.9p on A, and again with u as in (I3.35p . we infer 
from (I3.20p that x 1 —)■ (0|T“j,(x)|0) is in U°(r2) for all /i, G the 

regularized stress-energy VEV is continuous up to the boundary) fl^^l) . 

Now, let us £x a point x G hi and, making the assumptions fl3.8p . consider the 
functions u 1 —)■ (0|T“j,(x)|0); these are not only well defined, but even analytic on 
the half plane {u G C | 'Stu > d -|- 1} (Q). 

Let us recall that, according to Eq.s fl2.18p fl2.20p . the analytic continuation of 
(0|T“^(x)|0) at M = 0 determines the zeta renormalized VEV of the stress-energy 


^°This is slightly abusive, since staticity occurs only after taking the VEV; the alternative 
notation (0|T'“,,|0)(x) is more precise, but graphically disturbing and will not be employed in the 
sequel. 

^^Generalizing the previous considerations, one proves the following for any £ G N: with the 
assumptions (13.81) (resp., (Ih.Qp l. if > c?-|-£-|-1 the function x (0|T“j,(x)|0) is in C^(fl) (resp., 
in C^{n)). 

^^This fact can be established proving analyticity with respect to u of (x, x), of (x, x) 
and of the second order derivatives of at (x, x). We do not go into the details of the proofs, 
that will be given elsewhere m using rather obvious analyticity results for the operator functions 
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tensor; of course, the latter does not depend on t as well and we will write 


{0|f^j,(x)|0)^en = (0|f^;,(t,x)|0)^en 


(3.36) 


Of course, due to Eq.s (I3.311l3.33ll . the renormalized stress-energy VEV is determined 
by the “renormalized” functions 



RP 


u=0 


k“Z1;^(x, y) 


(3.37) 


■= RP 

2 u=0 V 2 


(3.38) 


(with z,w any two spatial variables), to be evaluated along the diagonal y = x. 
More precisely, we have 


(0|Too(x)|0)^en = 


(x, y) + ( 7 ) (d^^dye + V(x))Dl'7 (x, y) 


(0/ 


, (3.39) 


y=x 


(0|foj(x)|0)^en = (0|fjo(x)|0)^en = 0 , 
(0|Ty(x)|0)^en = (0|Tji(x)|0)^en = 



(3.40) 


(3.41) 


Let us remark that, if L*;!i±i(x, y) and d;z^Dv±i(x,y) have analytic continuations 
regular at n = 0, indicated hereafter with Zl^i(x, y) and dziuRi(x,y), one has 


^±i(x,y) = L>±i(x,y) , 

in 2 2 

dz^,D[''\x,y) = a^^L)i(x,y) 
2 2 


(3.42) 


for any choice of the mass scale k; clearly, in this case the renormalized stress- 
energy VEV is independent of k. On the contrary, an explicit dependence on k 
appears if the analytic continuations of Zlii±i(x, x) or (x, x) (or both) have 

a singularity at u = 0; this will occur in some specihc examples, to be considered in 
the subsequent Parts II and III. 

A connection between a regularized stress-energy VEV and a Dirichlet-like kernel 
is mentioned by Cognola, Vanzo and Zerbini [21] in a slightly different framework, 
where the held theory becomes Euclidean after Wick rotation of the time coordi¬ 
nate, and our operator A (in the spatial variables x) is replaced by the (spacetime) 
differential operator —du — A + V. 
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A variant of the previous results about the Dirichlet kernel and the regularized VEV 
(0|T“^(x)|0) can be formulated in the case of a slab. In this case hi = hli x 
with rii a domain in and di + d 2 = d; moreover, the potential V depends only 
on the coordinates xi G hli- In this situation we can express (0|T“j,(x)|0) in terms 
of the Dirichlet kernel associated to the operator Ai := —Ai + D(xi) (we defer to 
subsection 13.181 a comprehensive discussion of slabs configurations). 

In the following subsections we return to the case where D is an arbitrary domain 
in R'^ and we connect the Dirichlet kernel Dg to other integral kernels, in order to 
shed light on the analytic continuation of Dg (of course, these connections will also 
be useful, in their di-dimensional formulation, in the case of a slab). 

3.9 The heat kernel, the cylinder kernel and some variations. Let us con¬ 
sider again a strictly positive selfadjoint operator A in L‘^{Q) (that will be —A -|- V 
in the subsequent applications). For all t G [0, -|-oo), using the standard functional 
calculus we can dehne the operators 

e -^, . (3.43) 


These fulhll the following conditions: 


Jt 


moreover, due to the strict positivity of A , e~^ and are expected to vanish 

for t —)■ -I-CX3, in some sense that can be made more precise in terms of integral 
kernels. Let us now pass to the kernels 


= 0 , 


i=0 ^ ’ 

(3.44) 

-tv^ ^ Q ^ 

e-iVA 

= 1 ; 

(3.45) 


/L(t;x,y):=e-‘^(x,y) , T(t; x, y) = e-‘^(x, y) , 


(3.46) 


which can be expressed as follows in terms of a complete orthonormal set {Fk)k^ic 
of eigenfunctions of A and of the corresponding eigenvalues {coDkeic- 


A'(t;x,y)= / dke Ffc(x)Ffc(y) , 

Jk 

T(l;x,y)= f d/c Ffe(x)Ffc(y) . 


IK 


We note that 


(dt + ^x)^(t;x,y) = 0 , A:(0;x,y) = d(x-y) ; 

{dii- AA>T{i]^,y) = 0 , T(0;x,y) = d(x-y) . 


(3.47) 

(3.48) 

(3.49) 

(3.50) 
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In the above indicates the operator A acting on (t; x, y) and T(t; x, y) as func¬ 
tions of the X variable; Eq.s fl3.49p (13.501) follow, respectively, from Eq.s fl3.44p (13.451) . 
Besides, under minimal supplementary conditions one can prove that i^(t; x, y) and 
r(t;x, y) vanish exponentially for fixed x, y G hi and t —>■ -fcxo; we shall return on 
this in subsection 13.121 

In case (13.81) (i.e., A = —A -|- V with V smooth) we have di -|- ^x = — Ax -|- V (x) 

and dii — ^x = so Eq. (|3.49p contains a heat equation and Eq. 

(I3.50p a (d -I- l)-dimensionaI Laplace equation (with an external potential); note as 
well that iL(t;x, y) and T(t;x, y) fulhll the boundary conditions in the definition 
of A for X or y in dQ. For obvious reasons, K is called the heat kernel of A (even 
in cases where A is not of the form — A-l-I/); T is called by Fulling |1H] the cylinder 
kernel of A. It should be mentioned that, under the assumptions (13.8p . there are 
rigorous proofs that 

(t,x,y) !-)> iL(t,x,y), (t,x,y) H- T(t,x,y) are in C°°((0,-Fcx)) x x fi); (3.51) 

with the stronger assumptions (13.9p there is a result similar to (I3.5ip . with hi replaced 
by n. (These results can be proved by a slight generalization of Thm. 5.2.1 in [28] . 
For some related statements, see also Appendix |B] and [4T].I 

Again with the assumptions (I3.9p . A has pure point spectrum and possesses a com¬ 
plete orthonormal set of proper eigenfunctions (Ea,)^^! 2 , 3 ,..., for which the expansions 
(I3.47P (I3.48P hold with j^dk = One can prove that, for each t > 0, these 

expansions converge absolutely and uniformly on hi x hi with their derivatives of any 
order (see again Appendix El especially Eq. (IB.351) ). 

With the more general assumptions (13. 8 p we have A^ = A = A, which in turn 
implies similar relations for e~^, due to Eq. (13.6p . this gives 

Er(t;y,x) = A:(t;x,y) , T(t;y,x) = T(t;x,y) . (3.52) 

Needless to say, the heat kernel has been the object of intensive and detailed studies, 
even in a much more general framework than the one considered in the present paper; 
exhaustive analyses have been given, for example, by Berline et ah nn, Calin et al. 
[T8] . Chavel [22], Davies [28], Gilkey [53] and Grigor’yan [5l]. On the contrary, the 
cylinder kernel is a less popular object; it has mainly been investigated by Fulling 
and co-authors [m HHl ED]. Some considerations of Fulling (see, e.g., ISH) also 
involve the operator \fA associated kernel 

E(t;x,y) := (VA ^e"‘'^)(x,y) = f ^ Ffc(x)Ffc(y) , (3.53) 

Jk 

which we will refer to as the modified cylinder kernel, for reasons which become 
apparent hereafter (see Eq. (I3.54p j. Let us observe that the trivial relation = 
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can be reformulated in terms of integral kernels as 


r(l;x.y) = -a,r(t;x,y) ; 


(3.54) 


conversely, T can be determined as the primitive of —T which vanishes for t +cx3, 
that is 

/•H-OO 


r(t;x,y) = / dt' T(t';x,y) . 


(3.55) 


In some cases T is easier to compute than T, and some identities relating the cylin¬ 
der kernel T to the Dirichlet kernel Dg can be applied more efficiently if they are 
rephrased in terms of T (this situation will be exemplihed in the case of a wedge 
domain, to be discussed in the subsequent Part II; see Section 5 therein). 

With the assumptions fl3.8p or fl3.9l) . one could give for T some regularity results 
very similar to the ones illustrated previously for T (e.g.: statement fl3.5ip holds as 
well for T under the conditions (13.81) h 


Before moving on, let us consider the heat and cylinder traces; these are respectively 
dehned, for t G (0, +oo), as 


K(t) := Tre"*-^ , T(t) := Tre'^^ . (3.56) 


Assuming the above traces to exist, the general identity (13.4p for the trace of an 
operator B (here applied with either 13 = e~^ or B = e~^^) yields respectively 


K{t) = [ dx iP(t; X, x) , T(t) = f dxT(t;x,x). (3.57) 

Jn Jn 

In particular, in the case (13.9p where the eigenvalues of A are labelled by /C = 
{1, 2, 3,...} and the Weyl estimate (I3.inp holds, we have: 

H-co H-co 

K{t) = < +00 , T(t) = < +00 for all t > 0 . (3.58) 

k=l k=l 

Clearly enough, an analogous discussion could be made for the space integral of the 
diagonal, modihed cylinder kernel T(t;x, x); we omit this discussion for brevity. 

3.10 The case of a non-negative A. The heat and cylinder kernels. Let us 

remark that the heat and cylinder kernels can both be dehned even in case A is non¬ 
negative, without requiring strict positivity; by this we mean that <j{A) C [0, -|-oo), 
and that we are not assuming <j{A) C [e^, -f-oo) for any e: > 0. The non-negativity 
of A is equivalent to the existence of a complete orthonormal system {Fk)k£ic of 
(either proper or improper) eigenfunctions with corresponding non-negative eigen¬ 
values {uDkeic ^ 0). Most of the considerations of the previous subsection still 
hold, in particular Eq.s (I3.47P (I3.48p . 
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For example, if ^ = —A and 17 = R*^, then A is non-negative with eigenfnnctions 
Fk(x) = and eigenvalnes ui^ = |kp, labelled by k G R'^; the measnre 

dk on the set of labels is the usnal Lebesgue measnre of R'^. The eigenfnnction 
expansion fl3.47p of the heat kernel yields in this case the familiar resnlt 

; ( 3 . 59 ) 

moreover, the expansion fl3.48l) of the cylinder kernels gives the resnlt 


r(t;x,y) 


r( 


(j+l 

2 


)t 


d+l 
TT 2 


(F + |x-y|2) 


d+1 ’) 
2 


(3.60) 


which is a bit less popnlar and appears, e.g., in [49] . 

In some snbcases with non-negative spectrnm we can speak as well of the modihed 
cylinder kernel T(t;x,y). In fact, if 0 has zero spectral measnre (a fact holding 
when 0 belongs to the continnons spectrnm, bnt not holding when 0 is a proper 
eigenvalne), vCT can be dehned throngh the standard fnnctional calcnlns 

for selfadjoint operators. With minimal snpplementary assnmptions of regnlarity, 
the kernel T(t;x, y) := (dxlV^ e~^'^5y) makes sense, as well as its expansion 
fl3.53p in terms of a complete orthonormal set of (generalized) eigenfnnctions Fk 
with eigenvalnes (/cG/C); note that the assnmption of zero spectral measnre for 
0 is eqnivalent to the reqnirement that = 0 only on a zero-measnre snbset of K. 
In these sitnations we have again Eq. fl3.55p . describing the cylinder kernel T as the 
primitive of — T . 

For example, let us return to the case where A = —A and 17 = R'^, in which the 
spectrum (j{A) = [0, -|-cxd) is purely continuous. We have mentioned previously 
the eigenfunctions Fk and the eigenvalues (labelled by kGR'^), where Fk(x) = 
(27r)“'^/^e**^'^ and Uk = |k|; of course Uk = 0 only on a set of zero Lebesgue measure 
(consisting of the unique point k = 0). In this case, the expansion fl3.53p ( involving 
an integral in the Lebesgue measure dk) gives the following for d ^ 2: fl^^l ) 


^(t;x,y) 


r(^) 

27r^(F + |x - yP)^ 


(3.61) 


3.11 Connections between the cylinder kernel and a ((i-|-l)-dimensional 
Green fnnction. Due to the limited popularity of T it can be useful to connect 
this kernel to a more familiar object, namely a Green function, even though this 
requires to pass to d -|- 1 dimensions. All details of this construction are given in 
Appendix O where, as an example, this approach is used for a novel derivation of 
Eq. fl3.60l) not relying on the eigenfunction expansion fl3.48p . 

^^For d = 1 the right hand side of Ea. (l3.53p for T does not converge (not even distributionally); 
we take this as an indication that the modified cylinder kernel is ill-defined. 
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3.12 Behaviour of the heat and cylinder kernels for small and large t. 

The small t limit. The asymptotic expansion of the heat kernel of an operator 
A = —A + V on an open set C R'^ has been extensively studied (see, e.g, 
the work of Minakshisundaram and Pleijel m, or the already cited monographies 

on the heat kernel). From here to the end of this paragraph, we 
discuss the t —)■ O’*" behavior of the heat and cylinder kernels of .4, = —A+V under the 
assumptions fl3.9p (these could be generalized to get the same results, but we are not 
going to discuss this subject). As known from the previosly cited references, there is 
a unique sequence of real functions a„ : x —)■ R (n = 1, 2, 3....), usually referred 
to as HMDS (Hadamard-Minakshisundaram-DeWitt-Seeley) coefficients, such that 
for any N G {1, 2, 3,...} one has 


^(t;x,y) = 


(dvrt)'^/" 


|x-y|" 

e 4t 


N 


l + 5^a„(x,y)r + 0(t'^+i) 


n=l 


fort^0+. (3.62) 


In the above equation notice the factor iFo(t;x, y) := e ' , which is just 

the heat kernel associated to —A on R*^. In case H = 0, we have a„ = 0 for all n; 
thus, iF(t:x, y) = iFo(t: x, y)[l + 0(t°°)l (where the last term indicates a remainder 
which is O(t^) for each N e {1, 2, 3,...}). 

Along the diagonal y = x (for any V) Eq. fl3.62l) reduces to 

N 

1 + t'" + C>(t^+^) for t ^ 0+ 


A(t;x,x) = 


(47|-t)d/2 


n=l 


(3.63) 


where an(x) is shorthand for a„(x, x). The small t analysis of the cylinder kernel is 
more involved; however. Fulling proved (see, e.g., m) that its asymptotic behaviour 
along the diagonal y = x is as follows: there exist functions Cn, fn '■ ^ R 
(n = 0,1, 2,...) such that, for any N G {0,1, 2,...}, 


T(t;x, x) 


1 


N N 

J2en{^)i^+ /n(x)rinl + 0(t'^+dnt) 

n=0 n=d+l 

n—d odd 


for I —)■ O’*". (3.64) 


As pointed out in m, some of the functions e„, /„ (but not all of them) can be 
expressed in terms of the diagonal HMDS coefficients x i-G- a„(x) mentioned before. 


Before proceeding, let us remark that the heat and cylinder traces iF(t),T(t) (see 
Eq.s fl3.56p fl3.57p j are well-known to admit small t expansions analogous to those in 
Eq.s fl3.63p fl3.64p : see once more the references cited above. In particular, assuming 
again D to be compact and V to be smooth and bounded below, for t —)■ 0+ there 


hold 


A(l) 


1 

(47rt)c*/2 


N 

Ho/(D)+5^Aa"/' + 0(t'^) , 

71=1 


(3.65) 
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(3.66) 


T(t) 


1 

F 


N N 

Y,Ent^+ Fnt^\nt + 0(F+Mnt) 

n=0 n=d-\-l 

n—d odd 


{Vol{fl) denotes the volume of the spatial domain O). Notice, in particular, that the 
expansion fl3.65p for K{t) involves half-integer powers of t, whereas in expansions 
fl3.62p (13.6311 for the local heat kernel ii"(t; x, y) only integer powers of t appear; 
besides, let us stress that the real coefficients An, En, Fn in Eq.s fl3.65p fl3.66p are 
not just the integrals over the spatial domain O of the functions a„(x), e„(x), /„(x) 
of Eq.s fl3.63p (13.6411 . since boundary contributions arise as well. 


The large t behavior. Let us move on and note that, as anticipated in subsection 13.91 
both the heat and the cylinder kernel (along with their traces) vanish exponentially 
for large t, under minimal regularity conditions. As an example, we will prove the 
exponential decay of the cylinder kernel with the following assumptions: 

a) ^ is a strictly positive, selfadjoint operator in possessing a generalized 

complete orthonormal set {Fk)keK. of eigenfunctions with eigenvalues ul {uk>e>0)] 
the Efc’s are continuous functions on ff. 

b) for all points x, y G and t > 0, one has 


f(t;x,y):= f dk e |Ffc(x)| |Ffc(y)| <-Fcx) . (3.67) 

JK 

In this case, starting from the eigenfunction expansion (I3.18P we hnd 

l^(t;x,y)| ^ f(t;x,y) (3.68) 

at all points x, y G 14. On the other hand, after hxing r > 0 we see that 

T(t;x, y) ^ x, y) for all t ^ r (3.69) 

(this follows from the dehnition (I3.67p of T, noting that the inequalities Uk ^ £ > 0 
and i—T ^ 0 imply ^ The relations (I3.68P and 

(I3.69P give the desired result of exponential decay for T, since at all points x, y G O 
they imply 

|T(t;x,y)| ^ e-"(‘-")f(r;x,y) for all t ^ r . (3.70) 

In the above argument, condition (I3.67P plays a crucial role. This condition is 
fulhlled, for example, if ^ = —A -|- V and the assumptions (13. 9 p are satished. In 
this case we have a uniform bound T(t;x, y) ^ T(t) < -foo for all x, y G O and 
t > 0 (see Appendix [Bl Eq. (IB.37P ). so that Eq. (I3.70p implies 

|T(t; X, y)| ^ T(r) for all x, y G 14 and t ^ r . (3-71) 


The exponential decay of the heat kernel and of the traces of both the heat and 
cylinder kernels can be derived by similar considerations. For alternative approaches, 
see [28] or [54] . 
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3.13 The Dirichlet kernel as the Mellin transform of the heat or cylinder 
kernel. The results reported in this subsection are well-known; they were derived 
by Dowker and Critchley [30], Hawking [58] and Wald [101] and later reconsidered 
by Moretti et ah (see [IHITT] and citations therein). 

The representations of Dg mentioned in the title are useful in view of the analytic 
continuation with respect to s; they can be derived starting from the well-known 
relation (see [80], p.l39, Eq.5.9.1) 

+00 

di for all z G (0, -f 00 ), s G C with JJs > 0 . (3.72) 



This identity, along with the eigenfunction expansions fl3.18p . fl3.47p and fl3.48p of 
the Dirichlet, heat and cylinder kernels, yields 


^.(x,y) = 

^s(x,y) = 


r‘+oo 


r(s) 

1 

fM 


dtf A'(t;x,y) , 


b+co 


dtt^ T(t;x,y) . 


(3.73) 

(3.74) 


For example, Eq. (13.731) is derived via the following chain of equalities: 


^s(x,y) = 
1 


dk 

,2s 


r(s) 


X 

' + CXD 

di 


Ffc(x)Ffc(y) = dk 


^+OD 


r(s). 


dA:e-‘"EF,(x)Ffc(y) = 


dtE-'e-‘"^'Ffc(x)Ffc(y) = 

1 ,+00 (3-75) 

dtE-ii7(t;x,y) . 


'0 


Ik 


rW. 


In the second passage above, we used Eq. (13.721) with z = ; in the third passage, 

the exchange in the order of integration is justifed with arguments similar to those 
in [731171]. The derivation of Eq. (I3.74p is similar; in this case one has to resort to 
Eq. fl3.72p with z = ojk and s replaced by 2s. 

Eq.s fl3.73p fl3.74p state that the Dirichlet kernel Dg can be represented as the Mellin 
transform of either the heat or the cylinder kernel; there are analogous relations for 
the derivatives of the Dirichlet kernel, involving the corresponding derivatives of the 
heat and cylinder kernels. 

In the sequel we are especially interested in the case in which the integral repre¬ 
sentations fl3.73p and fl3.74p hold pointwisely, at a given pair of points x, y G D 
(including the case y = x). This occurs under the following three conditions: 

a) Dg(x,y) is well dehned at the given points x, y; 

b) K(t;x,y) and T(t;x, y) are also well-dehned at these points, for all t > 0. 

c) the integrals in Eq.s fl3.73p fl3.74p converge. 

With the assumptions fl3.9p . a) b) c) hold at all points x, y G D if 3fJs > d/2. In fact: 
for > d/2, Dg is continuous on D x D (see Eq. fl3.20p h so there is no problem 
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with its pointwise evaluation; for each t > 0, the functions K{i] , ) and T(t; , ) 
are continuous (in fact C°°) on x and again their pointwise evaluation is not 
a problem; due to the exponential decay of it'(t;x, y) and T(t;x, y), the integrals 
fl3.73p fl3.74p have no convergence problems for large t; for any x, y G 14, due to the 
t —?■ O’*" asymptotic expansions 03.621) 03.641) . the same integrals are both convergent 
for t close to zero if 3f4s > (i/2 ( 


14 


The pointwise representations 03.73P 03.74|) can be used as a starting point to build 
the analytic continuation of the function s i—)■ y) to a larger domain than 

the one where they are initially granted to hold (e.g., larger than the half plane 
{s G C I 3f4s > d/2} of case 03.9p ): we return to this point in the next two subsections. 
One could make similar considerations for the derivatives of Dg. For example, Eq. 
(I3T3D gives 

1 


P + CXD 


a“a^T>,(x,y) = 


r(.) 


dtF-'a“a^iF(t;x,y) ; 


(3.76) 


for any pair of multi-indices o;, jS. If we make the assnmptions 03.91) on ^ = —A -|- V 
and admit that the t ^ O’*" asymptotic expansion 03.621) can be derived term by term, 
1^. 03.76P holds pointwisely at all x, y G 14 if 3f4s > d/2 + j/2, where j := |q;| -|- |/5| 

(B- 

To conclude this subsection let us notice thah setting y = x in Eq.s 03.73P 03.74p 
and integrating over the spatial domain 14 ([/j); the relations 03.22p for the trace 
Trand 03.56P 03.57P for the heat and cylinder traces iF(t), T(t) allow us to infer 


Tr^-'* 


b + OD 


r(s) 


dt F-i K{i) ; 


(3.77) 


TtA-^ 


b+co 


r(2s) 


di T(t) . 


(3.78) 


Eq.s 03.77P 03.781) can be used to continue analytically the function s i—)■ Tr (x4“®); 
the situation is similar to the one outlined previously for the local counterparts of 
these equations, and will be reconsidered in the next two subsections. 


^■^For example, from (I3.62p it is clear that, in the limit i —>■ O’*', one has K{i,x,y) = 0(1) for 
y 7 ^ X and K{t,K,x) = 0 (l/F/ 2 ). so the integral in (13.731) converges for all x, y if Jfs > (i/2. 

^^In fact: for Jfs > d/2 +7/2, Dg is on O x (see Eq. (13.201) ). so there is no problem with 
pointwise evaluation of its derivatives up to order j; for each i > 0, the function K{i; , ) is C°° 
on O X 12, and again the pointwise evaluation of its derivatives is not a problem; the derivatives 
of d“d^i7(t;x, y) are easily seen to decay exponentially, so the integral (13.761) has no convergence 
problem for large i; deriving term by term the asymptotic expansion (13.621) we see that, for t —>■ O'*', 
d^d^K{t-,x,y) = 0(1) for y 7 ^ X and x, x) = (p(i/F/2+[i/2]^ which implies that, for 

5is > d/2 + j/2, the integral (j3.76|) is in any case convergent. 

Assuming that the order of integration can be interchanged for s is a suitable complex domain. 


29 
























































3.14 Analytic continuation of Mellin transforms via integration by parts. 

In the first part of this subsection, the analytic continuation via integration by parts 
will be presented for an arbitrary Mellin transform; in the second part, we will 
connect this general construction to the representation of the Dirichlet kernel (resp., 
of Tr as the Mellin transform of either the heat or the cylinder kernel (resp., 
of their traces). 

Let T : (0, +cxo) ^ C be a function of the form 

.F(t) = iw(t) (3.79) 

for some p G C and some smooth function "H : [0, +cx)) —)■ C, vanishing exponentially 
for t +cxo; consider the Mellin transform of i.e., the function 


fm(cT) 


+ 00 

di r-^ 


^(t), 


(3.80) 


dehned for appropriate cr G C. Due to the hypotheses on T, the integral in Eq. 
fl3.80p converges only for a G C with JJa > JJp and gives an analytic function of a 
in this region. However, integrating by parts n times (for any n G {1, 2, 3,...}) and 
noting that the boundary terms vanish (for JJa > JJp), we obtain 


911(a) 


(-ly 


r*+oo 


(a-p)...(a-p+n-l) Jq 


dtr-^+”-i^(t). 


(3.81) 


In consequence of the features of the function "H, the above integral converges for 
3?a > — n; thus, Eq. fl3.8ip yields the analytic continuation of the Mellin 

transform 91t(a) to the region 


{a G C I 3?a > — n} (3.82) 

from which the zeros of the denoninator in fl3.8ip must be removed; this gives a 
meromorphic function with (possibly) simple poles at the above zeros, which are 
the points 

a G {p, p - 1, ..., p - n + 1} . (3.83) 

Moreover, since the above results hold for any given n G {1,2,3,...}, they actually 
allow to determine the analytic continuation of 91t(a) to the whole complex plane 
with simple poles at the points a G {p, p — 1, p — 2,...} . 

As mentioned before, the above results can be employed to obtain the sought-for 
analytic continuation of the Dirichlet kernel Dg (treating x, y G D as hxed parame¬ 
ters) starting from its representations fl3.73p fl3.74p in terms of the heat and cylinder 
kernel, respectively. 
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More precisely, consider the case in which the heat or the cylinder kernel is given 
by a smooth fnnction of t rapidly vanishing at infinity, divided by a power of t; by 
this we mean that 

^(t;x,y) = ;^ iy(t;x,y) or T(t;x,y) = J(t; x,y) , (3.84) 

where p,q E H, H, J : [0, +cxo) x hi x hi —)■ R, and it is assnmed that (for hxed 
X, y G n) the fnnction t G [0,+oo) i—)■ if(t;x, y) or J(t;x, y) is smooth and rapidly 
vanishing for t — )■ +cxo. In these cases the integrals in the right-hand sides of Eq.s 
fl3.73p and fl3.74p converge for > p and 3fJs > q/2, respectively. Let ns mention 
that, with the assnmptions fl3.9l) . the heat kernel oi A = —A -f 1/ is as in Eq. fl3.84p 
with p = d/2 (recall the asymptotics fl3.62p ): again with the assnmptions fl3.9l) . the 
cylinder kernel of A is as in (I3.84p with q = d, provided that no logarithmic terms 
appear in the asymptotic expansion (I3.64p . 

Under the previons assnmptions, 1^. fl3.8ip . along with Eq.s fl3.73p fl3.74|) . gives the 
following for any n G {1, 2, 3,...} fj^ : 


^.(x,y) = 


^.(x,y) = 


(-1)' 


^+oo 


T{s){s-p)...{s-p+n-l) Jq 


dtU-^+"-iarR(t;x,y) ; (3.85) 


(-1)' 


'•- 1-00 


, --/ dtt"'*-''+'^-^a"J(t;x,y) . (3.86) 

T{2s){2s-q)...{2s-q+n-l) Jo i K, .yj K ) 

Comments analogous to the ones below Eq. fl3.81l) can be done for the above repre¬ 
sentations. More in detail, on the one hand Eq. fl3.85p gives the analytic continuation 
of the Dirichlet kernel Dg in the region {s G C | 3fJs > p — n} to a meromorphic 
function with simple poles at s G {p, p — 1,..., p — n -|- 1} ; on the other hand, Eq. 
fl3.86p gives the analytic continuation of Dg to the region {s G C | > (g — n)/2}, 

with (possibly) simple poles at s G {g/2, (g — l)/2, ..., (g — n -|- l)/2} . 

Of course, relations analogous to fl3.85|) and fl3.86p hold as well for the spatial deriva¬ 
tives of the Dirichlet kernel. 


In conclusion, let us stress that similar results can be deduced for the trace Tr A~‘^ 
(see Eq. fl3.22p L giving its analytic continuation to wider regions in the complex 
plane. For example, assume the heat trace has the form (compare with the hrst 
relation in Eq. fl3.84p ) 

A'(t) = 4 HU) , (3.87) 

for some p G R and some smooth function H : [0, -|-cx)) —)■ R, rapidly vanishing 
for t —)■ -|-cx); then, starting with Eq. fl3.77p and using the relations fl3.79113.8ip . we 

obtain Eq. (13.851) one uses Eq.s (13.79113.8fp with J^(i) = iir(t;x, y), p = p, 'H(t) = 
i7(t;x, y)) and cr = s. To obtain Eq. (13.8611 one uses Eq.s (13.79113.8fp with J^(t) = T(t;x, y), 
p = q, 7f(t) = J(t;x, y) and a = 2s. 
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obtain the following, for n G {1, 2, 3, 


Tr^-" = 


(-1)' 


f + OO 


r(s)(s-p)...(s-p+n-l) 7o 


dii 


Jn Tj 
s-p+n-1 ^ ^ 


cit" 


(t) 


(3.88) 


The above relation gives the analytic continnation of Tr^“® to the region {s G 
C I 9?s > p — n} to a meromorphic fnnction with simple poles at s G {p, p — 
— n + 1}. A similar resnlt can be derived nsing the cylinder trace T(t). 

3.15 Analytic continuation of Mellin transforms via complex integration. 

Another way to obtain the analytic continnation of the Mellin transform of a given 
fnnction is available (assnming the latter to fnlhll snitable conditions). Consider 
again the framework of the previons snbsection; this time the idea is to re-express 
the integral in Eq. fl3.80p as an integral along a snitable path in the complex plane. 
To this pnrpose, hrst consider the following identity concerning Mellin transforms. 
Let t h{t) be a complex-valued function, analytic in a complex neighborhood of 
[0, -|-oo) and exponentially vanishing for JJt —?• -l-cxo in this neighborhood; then 

^+CO p 

/ di T-^h(l) =-^ / dt T-^h(t) for s G C\{1, 2, 3,...}, 3fJs > 0 , 

Jo 2? sm(7rs) 

(3.89) 

where denotes the Hankel contour, that is a simple path in the complex plane 
that starts in the upper half-plane near -foo, encircles the origin counterclockwise 
and returns to -l-cxo in the lower half-plane (see Fig. [1] below). 





m 


In the right-hand side of Eq. fl3.89p . the complex power is dehned making 
reference to Eq.s fl3.15p fl3.16p : in the left-hand side, since t G (0,-|-oo), we dehne 
according to the standard convention fl3.14p . See Appendix |D] for the derivation 
ofEq. ffM . 

Assume now J- is as in Eq. fl3.79p with "H a complex function, analytic in a neigh¬ 
borhood of [0, -|-cxd) and exponentially vanishing for 3fJt -|-cxo; then, considering 
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the Mellin transform Tl{a) in Eq. fl3.80p and nsing Eq. fl3.89p with s = a — p and 
h = 71, we obtain 






2i sin(7r((T —p)) 


di . 


(3.90) 


In principle, Eq. (I3.9np holds nnder certain conditions: to ensnre existence of 9Jl(cr) 
as dehned by Eq. fl3.80l) we mnst reqnire > 3fJp, and the denominator sin(7r((j—p)) 
mnst be nonzero. However, the integral in Eq. fl3.90p converges for any a G C, so 
this eqnation yields the analytic continnation of the Mellin transform J\A(a) to the 
whole complex plane, possibly with simple poles for 


a G {p, p- 1, p-2, ...} , (3.91) 

dne to the vanishing of the sine fnnction in the denominator (Q). 

Recall once more that, according to Eq.s fl3.73p fl3.74p . the Dirichlet kernel can be 
expressed as the Mellin transform of either the heat or the cylinder kernel; so, the 
above resnlts on the analytic continnation via contonr integration can be applied to 
Zls(x,y) (for hxed x,y G H). More precisely, snppose that either the heat or the 
cylinder kernel has the form fl3.84p : 

^(t;x,y) = ;^ i7(t;x,y) or T(t; x, y) = J(t; x, y) , 

with p, g G R and snitable fnnctions H, J : [0, +oo) x H x H R. Assnme these 
fnnctions to have extensions 77, J : 7/([0, +C)o)) x H x H —)■ C, where 7/([0, +C)o)) C C 
is an open neigbonrhood of the interval [0, +oo) and, for hxed x, y G H, the fnnction 
t G W([0, +cxd)) I—)■ i7(t;x,y) or J(t;x,y) is analytic and exponentially vanishing 
for —)■ +CXD. Making these hypoteses and nsing Eq. (I3.9np along with Eq.s (I3.73p 
(I3.74p . we obtain (pi ) 


^s(x,y) 


3-i7r(s-p) 


2ir(s) sin(7r(s-p)) 


[ dtE-i77(t;x,y) ; 


(3.92) 


^.(x,y) 


g-i7r(2s-(j) 

2ir(2s) sin(7r(2s —g)) 


dtR'*-'T(t;x,y) 


(3.93) 


^®By inspection of the denominator, it would seem that also the points a S {p+1, p+2, p+3, ...} 
are singular, but we know this is not the case since the original expression (13.801) for 9Jl(cr) is regular 
at these points. The reason for the apparent contradiction lies in the fact that the integral over 
the Hankel contour in Eq. (13.9011 vanishes for the above mentioned values of a (as can be easily 
checked via the residue theorem recalling the properties of J^), thus yielding an indeterminate form 
oo • 0. 

^®One proceeds as in Footnote [T71 using Eq. (13.901) in place of Eq. (13.811) . 
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Due to the comments after Eq. fl3.90p . both the above identities yield the analytic 
continuation of the Dirichlet kernel to a meromorphic function on the whole complex 
plane; more precisely, the analytic continuations obtained via Eq.s fl3.92p and fl3.93p 
may have simple poles respectively for s E {p, p — l, p — 2, ...} \ {0, —1, —2,...} and 
s G {g/2, (g — l)/2, (g — 2)/2, ...} \ {0, —1/2, —1, —3/2,...} (as readily understood 
analysing the denominators in the right-hand sides of the cited equations). 

In the subcases where p, g G Z = {0, ±1, ±2,using trivial trigonometric iden¬ 
tities and recalling that r(s)r(l — s) sin(7rs) = vr for any s G C (see [80], p.l38, 
Eq.5.5.3), Eq.s (13.921) fl3.93|) can be rephrased as 

D.{x,y) = ~ dt t-‘g(t; X. y) ; (3.94) 

B.(x,y) = iX,y) . (3,95) 

In these subcases the integrals along the Hankel contour can be computed straight¬ 
forwardly for integer and half-integer values of s, respectively, by means of the 
residue theorem; for example, for s = —n/2 and n G {0,1, 2,...}, Eq. fl3.95p yields 

(x, y) = (-1)" r(n+l) Res(t-("+') T(t; x, y); o) . (3.96) 

We can obtain a variant of Eq. fl3.95p . giving the analytic continuation of the 
Dirichlet kernel Dg in terms of the modihed cylinder kernel T (recall Eq.s fl3.53p 
fl3.55p h To this purpose, we assume T to admit a meromorphic extension in t to a 
neighborhood of [0, +cxd), having a pole in t = 0 and rapidly vanishing for 3fJt —)■ -|-cx); 
then, expressing the cylinder kernel T(t;x, y) in Eq. fl3.95p as —cltT'(t; x, y) and 
integrating by parts, we obtain 

B.(x, y) = - e-^‘"rp-2s) ^ f(t ; X, y) (3,97) 

(note that no boundary contribution arises, due to the rapid vanishing of T(t;x, y) 
for —)■ -|-cxd). Again, for half-integer values of s we can compute explicitly the 
analytic continuation fl3.97p by means of the residue theorem; to be more precise, 
for s = —n/2 and n G { — 1, 0,1, 2,...} we have 

D_5(x,y) = (-l)"+'r(n+2)Res(l-("+2)T(t;x,y);0) . (3.98) 

Relations similar to the ones obtained above hold for the spatial derivatives of the 
Dirichlet kernel, allowing in turn to determine their analytic continuations. 

To conclude, let us mention that similar results hold as well for the trace Tr^“®; 
these are obtained using the representations fl3.77p fl3.78p in terms of the heat and 
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cylinder trace K{t),T{t), and assuming suitable features for the latters. In partic¬ 
ular, if the map t i—)■ T(t) admits a meromorphic extension to a neighborhood of 
[0, -|-oo) which only has a pole singularity at t = 0 and vanishes exponentially for 
3?t —>■ -|-oo, for n G {0,1, 2,we have 

Tr^V2 = (-l)^r(n+l)Res(^r(’"+^)T(t);0) . (3.99) 

3.16 Other kernels, and their relations with Dg. In this section we are con¬ 
sidering a number of integral kernels connected with the zeta regularization of the 
stress-energy VEV; attention is mainly focused on the kernels used in the subsequent 
applications (including the subsequent Parts II,III and IV). However, it would be 
against the spirit of this section to ignore completely the resolvent kernel, i.e., the 
kernel of the operator {A — A)“^, where .4, is a selfadjoint operator and A G C is 
outside the spectrum a{A). Under appropriate conditions (in particular, the strict 
positivity of A), the powers A~^ can be related to suitable contour integrals in¬ 
volving the resolvent [95]; this fact can be restated in terms of a relation between 
the corresponding kernels. This possibility will not be considered here and in Parts 
II-IV, but we plan to recover it elsewhere. 

3.17 The case of product domains. Factorization of the heat kernel. Let 

us consider the case where A = —A -|- V and the spatial domain H C has the 
form 

U = ffi X ffs (3.100) 

with Ha, for a G {1, 2}, indicating an open subset of (di + d 2 = d); in this case, 
points of H will be written as 

x = (xi,x 2 ), y = (yi,y2) (3.101) 

etc., where x^, G (a G {1, 2}). In addition to Eq. fl3.100p . we assume that the 
external potential has the form 

I/(x) = I/i(xi) + U 2 (x 2 ) (3.102) 

and that the boundary conditions specihed on dfl = (dHi x H 2 ) U (Hi x 8 ^ 2 ) arise 
from suitable boundary conditions prescribed separately on dHi and 8 ^ 2 , in such a 
way that, for a G {1, 2}, the operator 

A:=-Aa + U(xa) (3.103) 

(Aa the Laplacian on H^) is selfadjoint in L^(Ha). Moreover, each one of these 
operators is assumed to be strictly positive or, at least, non-negative. 
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In the situation described above, the Hilbert space L^(f2) and the fundamental 
operator A := —A + H(x) acting therein can be represented, respectively, as 

(g) , A = Ai(S)1 + 1®A2. (3.104) 

Because of the assumptions we have made, each of the two operators Aa {a G 
{1,2}) possesses a complete orthonormal system of eigenfunctions {Fa^ka)ka&iCa with 
eigenvalues as for the fundamental operator A, we see that it has a complete 

orthonormal set of eigenfunctions of the form 

Ffc(x) := Fi,fci(xi)F 2 ,A: 2 (x 2 ) for k = {kiA 2 ) G /Ci x /C 2 (3.105) 

and that 

AFk = ujlFk , ujI = ^l^k^_F^ 2 ,k 2 ■ (3.106) 

Or course (j{A) = (j{Ai) + a{A 2 ), so that A is non-negative; besides, A is strictly 
positive if so is one at least between Ai and A 2 - 

In the product case under analysis, a number of interesting facts occurs for the 
integral kernels associated to A and Ai, A 2 - The most elementary of these facts is 
the factorization of the heat kernel; by this we mean that the kernels 

^(t;x,y) := (e"‘^)(x,y) , /^^(t; x„, y„) := (e"‘^“)(x«, y„) (aG{l,2}) (3.107) 

are related by 

A(t;x,y) = Ai(t;xi,yi) A 2 (t;xi,X 2 ) , (3.108) 

a fact that is made apparent by the eigenfunction expansion fl3.47p and by Eq.s 

(ITTO (ITT06D . 

In passing, let also mention that an analogous relation can be easily derived for 
the heat trace; writing K{t), Ka{t) for the heat traces of A and Aa (a G {1,2}), 
respectively, we obtain (Q 


K{t) = K,{t)K2{t) . (3.109) 

In the present subsection we have analysed the case of a product conhguration with 
two factors; as a straightforward generalization, one can consider a product with 
an arbitrary number of factors. Examples of such multiple products will appear in 
Parts III and IV. 

fact, Eq.s (13.1001) (13.1081) and the relations (I3.56P (I3.57P allow us to infer the following chain 
of equalities: 

K{i)= ( dxiy(t;x,x)= / dxi iyi(t;xi,xi) / dx 2 ^ 2 ( 1 ; X 2 , X 2 ) = iCdt) ^ 2 ( 1 ) . 

J ^ J J ^2 
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3.18 The case of a slab: reduction to a lower-dimensional problem. By 

definition, we have a slab if 

n = nixR‘^\ V = V(xi) (3.110) 

with fli a domain in (di + d 2 = d), and if the boundary conditions prescribed 
for the field refer to its behaviour on dfli x Clearly, a slab is a subcase of 

the general product case discussed in the previous subsection, with and 

1/2 = 0 . In this subcase the relevant operators are A = —A + V (xi) acting in 


Ai —Ai + V (xi) 


(3.111) 


acting in L^(f2i), and A 2 ■= —A 2 acting in 

The operator Ai has its own eigenfunctions Fi fcj(xi) = l^fci(xi) and eigenvalues 
uf = zul^ {ki G /Cl); we assume Ai to be strictly positive, so that zuk^ ^ £ 
for some e > 0. Of course, — A 2 is non-negative with eigenfunctions T 2 ,k 2 (x 2 ) = 
and eigenvalues ti; 2 k 2 ~ for k 2 G R*^^. 

In the sequel we write //^(xi, X 2 ; yi, y 2 ) for the Dirichlet kernel of A at the points 
X = (xi,X 2 ) and y = (yi,y 2 ); will be the Dirichlet kernel of Ai. 

On the one hand, according to the general results (13. 3 1113.331) . the regularized VEV 
(0|T“^(x)|0) is determined by Dirichlet kernel Ds(x.,y) and its derivatives evaluated 
on the diagonal y = x. The main intent of this subsection is to express Dg and its 
derivatives in terms of the reduced kernel at all points of the diagonal y = x 
(and, in fact, on an even larger domain). The starting point towards this goal is the 
identity 

^^(xi,X 2 ; yi,y 2 ) = yi; |x 2 - y 2 n , (3.112) 

involving a function Dg : Di X Di X [0, -t- 00 ) C, (xi, yi, q) i-G- D^(xi, yi, q). This 
function and its partial derivatives with respect to g, at g = 0, are completely 
determined by the kernel according to the following rules: 

T),(xi,yi;0)= yi) for s G C, > f ; (3.113) 


d^Dg 


(xi,yi;0) = 


(-l)’^r(s-f-n) 


A) 


D a, 

s-^-n 


(xi,yi) for sG C, 3fJs > ^-|-n. (3.114) 


dqu '' ^ (47r)‘^i/2 qn 

We defer to Appendix [E] the derivation of Eq.s fl3.112113.114p . In particular, for 
the derivatives involved in Eq.s (13.3 1113.35]) on {0|T“j,(x)|0) we obtain the following 
expressions (where, for simplicity of notation, we write (x, y) for (xi, X 2 ; yi, y 2 )): 


T>2i±i(x,y) 


r( 


R—^2^1 '' 

2 2 


(4vr)<i2/2r( 


uzLl '' 
2 ^ 


D 


( 1 ) 


(xi,yi) 


(3.115) 


y=x 
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yi=xi 














d i iDu+i(x,y) = d i iDu+i(x,y) = d i jDu+iiyi.y) 


= 0 


y=x 


for (a, 6) = (1, 2) or (a, h) = (2,1) and i G {1, ...da}, j e {1,..., 4}; 


r( 


U—d2 + l ^ 


for z,w & {x, ?/} and z, j G {1,..., di} ; 


yi=xi 


5^i^id}i+i(x,y) = - d i jDu±i{x,y) = - 9 i ,£>.+1 (x,y) 

^2^2 2 y=x ^2^2 2 y=X ^2^2 2 


Y ^ u—d2—l ^ 


(47r)<^2/2 2r(2i±l) 


y=x 

for ij G {1, ...,^ 2 } • 


(3.116) 


(3.117) 


(3.118) 


yi=xi 


Relations (13.115113.1181) are derived assnming JRm > ^2 + 1, bnt it follows from them 
that the analytic continnations in u of the Dirichlet kernel, of its rednced analogue 
and of their derivatives fulhll the very same relations. 

Let us remark that the left-hand sides of the above equations depend in principle 
on X = (xi, X2), while the right-hand sides only contain Xi; this confirms the expec¬ 
tation that the stress-energy VEV ought to be independent of the variable X 2 , due 
to the symmetry of the slab conhguration under translations regarding this variable 
alone. Finally, using Eq.s fl3.115113.118p and fl3.3HI3.33|l . it can be easily checked 
that the components of the regularized stress-energy tensor VEV also fulhll 


(0|7}}(x)|0) = 0 for z,j G {di + 1,..., d}, j ] 
(0|7}“(x)|0) = (0|7;“(x)|0) = 0 for zG{l,...,di}, jG{di + l,...,d} . 


(3.119) 


4 Total energy and forces on the boundary 

We refer again to the general framework of Section O where a quantized scalar held 
on a spatial domain hi and the associated stress-energy tensor VEV are considered; 
we recall that A indicates the fundamental operator —A -|- R(x) acting in L‘^{VL). 
From now on, we indicate with da the area element on the boundary dVL, and use 
for n the standard Lebesgue measure dx. We also write n(x) = (n^(x))£=i_,,,_d for 
the outer unit normal at a point x G 912; this is assumed to exist everywhere (which 
happens if 912 is globally smooth) or almost everywhere (which happens if 912 has 
edges or corners). 

4.1 The total energy. As anticipated in subsection l2.6L the zeta-regularized total 
energy can be dehned as 


[ dx(0fe(x)|0) , 
Jn 


(4.1) 
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provided that the above integral converges for m in a suitable complex domain; using 
Eq. (13.3111 . the above definition yields 


= + I dxD^{x,y) 


+ 


y=x 


+ k" --p/ <ix (a“s„<+r(x))r>s+.(x,y) 


(4.2) 


y=x 


On the other hand, the eigenfunction expansion (13.1811 for the Dirichlet kernel (here 
2 


used with s = gives 


dx 


In 


{d‘‘dyt + r (x)^ Da+1 (x, y) 


(4.3) 


y=x 


dk 


I I dx (+ l/(x)Ffc(x)Ffc(x)) = 

IK do ^ ^ 

dx (Ffc(x)(-a^(9£ + l/(x))Ffc(x))+ f da(x)Ffc(x)n^(x)(9£Ffc(x) 

I an 

fill 


r dk 

I ! 

where, in the last step, we have integrated by parts 
To go on we note that {— d^di + V)Fk = AFk = ojlFk and Ffc(x)n^(x)c?£Ffc(x) = 
Ffc(x)^^(x) = , where we have introduced the normal derivative 


y=x 


d/dn := n^d^e . Substituting into Eq. (14.3p and summing over fc G /C, we obtain 


dx 


In 


dyi + V(-x))Du+i(-x.\) 

d 


y=x 


(4.4) 


y=x 


= / dxT)«^(x, y) + / da{x.) —— Du+i (x. v) 

Jn ^ y=x Jan i^^y ^ 

inserting this result into Eq. (14.2p . we conclude 

+ B^ , 

where we have introduced the regularized bulk and boundary energies 


:= - / dxD^(x,x) = - Tr.A 


1 —Ti, 
2 


(4.5) 


(4.6) 




Here and in similar situations, whenever we speak of an integration by parts we refer to the 
identity 


>n Jan 

holding for all sufficiently smooth functions /, g . 


[ dyi{dif)g= [ dafgnt- f 
Jn Jan Jn 
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5 “ := K' 


y=x 


(4.7) 



d 

dUy 


^sM(x,y) 


The derivation of the above result is a bit formal since, in general, one cannot grant 
convergence of the integrals dehning and i?“, for suitable values m G C . 

As for the bulk energy, it is easy to give an example in which hniteness is granted for 
appropriate u. To this purpose let us consider the case fl3.9|) (involving a bounded 
domain with Dirichlet boundary conditions); in this case, recalling the Weyl esti¬ 
mates (13.101) . we conclude 


is hnite if > d-|-l . (4.8) 

As for the regularized boundary energy B'^ let us mention that, for bounded. 


= 0 under Dirichlet or Neummann boundary conditions on dVl (4.9) 


(since in the Dirichlet case we have Ds(x, y) = 0 for xg 9D and all y, while in the 
Neumann case g^Ds(x, y) = 0 for y G 9D and all x). 

Applying the above considerations in the case of an unbounded domain requires 
much caution. On the one hand, E'^ can be inhnite for all m G C; on the other 
hand, in the dehnition (14.71) of it might be necessary to intend the integral 
fg^ da as lim^^+oo Jg^^ da, where is a sequence of bounded subdomains 

such that VLi C (for any I G {0,1, 2,..}) and = D (note that this limit 

could either be inhnite or even fail to exist). 

If 77“ exists hnite for u belonging to a suitable open subset of C and it is an analytic 
function of u on this domain, a renormalization by analytic continuation can be 
implemented; in general, following the extended version of the zeta approach, we 
dehne the renormalized bulk energy as 


Even j^p 


E^ = RP 


u=0 


n ^ . l-rx 

, — Tr 2 

.=oV 2 


(4.10) 


P_ 

When the analytic continuation of Tryl^“ 
scription is reduced to 


^ren _ 


u=0 


is regular up to u 


Tr ^^/2 


0 the above pre- 


(4.11) 


(of course Tryl^/^ indicates the analytic continuation of Tryl^z" at m = 0). 

In a similar way one can dehne the renormalized boundary and total energies as 


pren pp 


gren pp 


5“ ; 

(4.12) 

. 

(4.13) 


4i=0 
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An alternative definition of the renormalized total energy could be 

[ dx(0|foo(x)|0),en. (4.14) 

Jn 

This possibility, which is considered rarely in this series of papers, is not granted to 
be equivalent to (Id.ldh : for example, it may happen that the integral in the right- 
hand side of Eq. fl4.14p diverges, while the prescription fl4.13p always gives a finite 
result by construction. For a comparison between the alternatives fl4.13p fl4.14p . see 
the final lines of subsection 16.41 (dealing with a field on a segment, for several types 
of boundary conditions). 

4.1.1 Reduced energy for a slab configuration. Let us consider the slab 
configuration introduced in subsection I3.18P so that 12 := ffi x potential V 

depends only on xi G 12i, and the boundary conditions regard only dl2i x 122. We 
already observed in the mentioned subsection that the regularized VEV (0|T^i,|0) 
depends only on Xi G 12i (and not on X 2 G R'^^); so, the integral in Eq. fl4.ip defining 
the total energy diverges due to an infinite volume factor. 

As a matter of fact, when dealing with a slab configuration one usually considers in 
place of the total energy the reduced total energy this is the total energy per 
unit volume in the “free” dimensions, i.e.. 


er := [ dxi (0|f;„|0) . (4.15) 

Recalling Eq.s fl3.115113.llSp and using some well-known identities regarding the 
gamma function, we infer 


= 


(47r)'^2/2r(^ 


^ u-l + d2 
y4:{u-l-d2) 


dxiZ2i^2d,_i(xi,yi) 


+ 


yi=xi 



) [ dxi 

(d^‘dye+V{:)Ci))D^^l^^+, (xi, yi) 

\ 


' Jfh 

_ '2 

yi=xi J 


(4.16) 


Concerning the second term above, we can express the reduced Dirichlet kernel 
in terms of the eigenfunctions and the eigenvalues (rufcJfcjgACi 

2 

of the reduced operator Ai = — Ai -|- R(xi) and integrate by parts as in the general 
setting; working as in the derivation of Eq.s (I4.3p (14.4p and keeping in mind that 
AdSki = we obtain 


dxi 


’ Qi 


dye + V{^i))D (xi, y 1 ) 


yi=xi 


JUl 2 


+ 


d 


yi=xi 


li —(^2 + 1 (xi,yi) 

IdQi 2 


(4.17) 


yi=xi 
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In conclusion, we have a result similar to Eq. fl4.5p : 


8 


U 

1 


+ B\ 


(4.18) 


where we have introduced the regularized reduced bulk and boundary energies 


: = 


^u-p^ n-d2-l 

2 (dvr 


p/ u-a2-i \ p 


r( 


u—d2 — l'' 


2 (47r)'^2/2 


Tr^ 


^2 + 1—1/ 


(4.19) 


:= 


K' 


r( 


U — d2 + l'' 


u+1'' 






d 


/ da(xi) 

Ian, dn- 


D 


yi 


( 1 ) 

li —<^2 + 1 
2 


(xi,yi) 


(4.20) 


yi=xi 


The considerations of the previous subsection about convergence of the bulk and 
boundary energies have obvious analogues for the reduced energies i?“. 

Of course, the reduced bulk and boundary energies are renormalized in terms of the 
analytic continuation (or, possibly, of its regular part) at m = 0. 


4.2 Pressure on the boundary. In this subsection we are interested in the pres¬ 
sure p(x) = (pi(x))j=i^...^d, i.e., the force per unit area produced by the field inside 
O at a point x on the boundary dVt. A possible characterization is the following: 
we first introduce, for Iblu large, the regularized pressure p“(x) of components 


P“(x) := (0|^“(x)|0)n^(x) for z G {1,..., d} ; (4.21) 


then, we define the renormalized pressure at x setting 


pr(x) := RP 


u=0 


.pr(x) 


(4.22) 


where RP\u=o indicates the regular of the analytic continuation evaluated at m = 0 
(of course, if the mentioned continuation is regular up to m = 0, the above pre¬ 
scription reduces to p[®"'(x) := pf{x.)\u=o, meaning that the analytic continuation at 
u = 0 has to be considered). 

It is important to point out that this is not the only reasonable definition for the 
renormalized pressure at a point x G ; another possibility is 


Pr”(x) := ( , lim (0|fij(x')|0)^en) ri^(x) . (4.23) 

In few words: in the approach (14.21114.2^ . one stays at a point on the boundary, and 
performs therein the renormalization; in the approach fl4.23p . one renormalizes at 
points inside hi, and then moves towards the boundary. Notice that both approaches 
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require the existence of the normal n(x) (and thus lose meaning on edges and corner 
points of Sfl). 

As a matter of fact, the prescriptions (14.21114.2^ and fl4.23p do not always agree. 
The approach (14.221) (possibly, in the restriced version) gives by construction a hnite 
pressure; on the contrary, this is not granted for the alternative prescription (I4.23p . 
As an example, in Part II of this series of papers we discuss the case where the spatial 
domain is a wedge; in this case at all boundary points not in the edge, where the 
normal is well dehned, the pressure dehned according to Eq. (14.231) diverges. We 
conjecture that, in general, at points x G dQ where the normal is well dehned and 
the approach (I4.23p gives a hnite pressure, the result obtained according to the latter 
prescription agrees with the renormalized pressure dehned by Eq. (14.221) : in fact, 
this happens in all the examples analysed in this series of papers. 

In the rest of the present section, our analysis of the boundary forces will mainly 
refer to the approach (14.21114.221) . 

In applications, one often considers a situation where a quantized held is present 
both inside and in the complementary region := R'^ \ . In this setting the 

force per unit area acting on the boundary is the resultant of the pressure produced 
by the held inside fl, on the one hand, and by the held inside on the other; the 
renormalized versions of both these observables can be computed, separately, using 
either one of the two approaches mentioned before. 

4.3 Explicit expression for the (regularized) pressure. Let us stick to the 
viewpoint (14.21114.22)) : in order to implement it, we use Eq. (I3.33P for the regularized 
stress-energy tensor that gives the following, for x G dQ (and n(x) well dehned): 


PiW = 


I -^)^p (Ds^(x,y) - (d*'dj^^ + I/(x))L)«M(x,y)) + 


~^ ( ( 2 ^^ ) D (x, y) 


(4.24) 


n-^(x) 


y=x 


To go on, let us restrict the attention to the case of Dirichlet boundary conditions; 
then, only the terms involving mixed derivatives (both with respect to x and y) of the 
Dirichlet kernel yield non-vanishing contributions on the boundary dQ. Moreover, 
the terms proportional to ^ in Eq. (I4.24p can be shown to vanish, so that 


p-{x) = 


- ^ 6ij ^ d^iyj ) (x, y) I nh x) 


1 

2 


(4.25) 


y=x 


see Appendix [F] for the proof. As a hnal step, the analytic continuation of at 
M = 0 must be considered. 


4.4 Au equivalent characterizatiou of boundary forces. In the literature, 
forces on dfl are often characterized by a diherent approach, which does not require 
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the knowledge of the full stress-energy tensor; see, e.g., the monographies by Bordag 
et ah Milton [72] and Plunien et ah [82]. In this approach one considers 

a variation of the spatial domain controlled by a real parameter, and dehnes the 
pressure in terms of the derivative of the bulk energy with respect to the mentioned 
parameter. For example, if is a parallelepiped (0, a) x (0, b) x (0, c) one could 
consider the variation of the length of any one of its sides, say a; it is customary to 
dehne the total force on the face = a} as the derivative of the bulk energy with 
respect to a, with the sign changed. 

The idea that boundary forces are related to the variation of has been typically 
presented in simple examples like the previous one; it can be of interest to propose 
a general formulation of this idea, and to compare it with the characterization of 
boundary forces given in subsection 14.21 via the stress-energy tensor. 

For the sake of dehniteness, let us consider the case where is a bounded domain in 
and Dirichlet boundary conditions are prescribed on dfl] besides, let & : —?■ 

R'^ be a vector held. We assume fl, its boundary dfl and © are regular enough to 
permit the subsequent calculations. 

First of all, consider the family of diffeomorphism 

5, ;R‘^^R'^, 5,(x) :=x + e6(x) , (4.26) 


labelled by a small parameter e > 0, that will be ultimately sent to zero. For any e, 
the spatial domain 

a := S,{n) (C R'^) (4.27) 

can be regarded as a deformation of the initial domain 12. 

Of course, a relation analogous to fl4.6|) holds for the regularized bulk energy 
associated to the deformed domain 12^, i.e.. 


E 


U 

e 


2 


fce/c 


1 —n, 
2 


(4.28) 


in the above (a;^^)fcgyc denote the eigenvalues of the fundamental operator As, that 
is the operator —A + V acting on the Hilbert space L‘^{Qs) (with Dirichlet boundary 
conditions on dfls) instead of L^(12). 

Let us now consider the expansion of 27“ to the hrst order in e, which describes the 
variation of the regularized bulk energy under the deformation 04.261) 04.271) of the 
space domain. Due to Eq. 04.28p . the calculation of this expansion can be reduced 
to the hrst order expansion of the eigenvalues Us,k', this can be done by standard 
perturbation techniques, as well known from the classic work of Rellich [85]. As 
illustrated in Appendix [F] the conclusion of this analysis is 


E^ = E^-e{l-u) [ da(x) ©*(x) p“(x) + 0{e^) (4.29) 

Jan 
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where = (p“) is the regularized pressure, given by Eq. fIF.lip . Eq. fl4.29p can be 
used for an alternative, but equivalent dehnition of the regularised pressure; from 
this viewpoint, the regularized pressure held is the unique vector held on dQ such 
that fl4.29p holds for each one-parameter deformation of the form fl4.261l4.27p for the 
domain 

Let us now perform the analytic continuation up to u = 0, assuming that no pole oc¬ 


curs at this point; E'^ 


u=0 


and 


u=0 


and Eq. fl4.29p yields the relation 


are the renormalized bulk energy and pressure. 


^ren ^ ^ren _ ^ ©*(x)p["^(x) + O(e^) . ( 4 . 30 ) 

Jon 

This is the result anticipated at the beginning of this subsection: a characterization 
of boundary forces in terms of the of the bulk energy variation under deformations 
of the domain. We already mentioned the frequent use of this idea in the literature, 
for particular choices of 

4.5 Integrated forces on the boundary. Let us now discuss the evaluation 
of the integrated force '^o acting on an arbitrary subset D of the spatial boundary 
(D C dfl; possibly, D = dfl). As in the case of the pressure considered in subsection 
14.21 we can give several alternative dehnitions of this quantity. As a hrst possibility, 
we introduce the regularized total force on D (for large 3fJM) 

5“ ;= / da(x) p“(x) , (4.31) 

Jo 

where p“ = (p“(x)) indicates the regularized pressure (see Eq. fl4.2ip i: then, we 
dehne the renormalized total force on D as 


do"" ■■= RP 


w=0 




(4.32) 


(clearly, when there is no pole in u = 0, the above prescription reduces to 
5ol«=0) meaning that the analytic continuation in u = 0 has to be considered). 
Another possibility is to put 


ro" := [ da{^) p^^"(x) , (4.33) 

Jo 

where p^^^Px) = (pf®"'(x)) is the renormalized pressure, dehned according either to 
Eq. fOTll or to Eq. fITTHll . 

Similarly to what we pointed out in subsection 14.21 for the pressure, in general the 
two alternatives fl4.32p fl4.33p give different results; in fact, the prescription fl4.32p 
always gives a hnite result for 5^"', while fl4.33p can give an inhnite result. 
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In conclusion, let us stress that for the integrated force there hold comments anal¬ 
ogous to the ones at the end of subsection 14.21 when a quantized held is present 
both inside and in the complementary region := R'^ \ . In this case the total 

force on any subset D C dQ is given by the resultant of the forces corresponding, 
respectively, to the held inside and outside the spatial domain . 

4.6 A comment on some previous “anomalies”. In the previous subsections 
we have pointed out that the renormalized versions of the total energy, of the pressure 
and of the integrated forces on the boundary can be dehned according to diherent 
prescriptions, which in general are not equivalent (see Eq.s fl4.13p fl4.14p . fl4.22p 
fl4.23p . (14.3211 (14.3311 and the considerations in the corresponding subsections). 

In consequence of this, there arise unavoidable ambiguities, or anomalies, when 
talking about the renormalized observables mentioned above. For example, we have 
mentioned previously the possible non-equivalence of the alternatives (I4.13P (I4.14p 
for the total energy and (I4.22p (14.2311 for the pressure recall that it may 
happen that the prescriptions (I4.14p and (I4.23P give inhnite results for and 
due to boundary singularities of the stress-energy VEV which make divergent the 
integral dx (0|Too(x)|0)ren or the limit limx/eo,x'^x(0|Tij(x')|0)renn^(x) (x G dfl, 
i e {l,...,d}). 

On the other hand, such boundary singularities of the renormalized stress-energy 
VEV are not a specihc consequence of the zeta regularization; indeed, they also 
appear if one uses point-splitting, as indicated by the very systematic analysis of 
Deutsch and Candelas [29]. For the moment, the above mentioned anomalies must 
be accepted as a problematic aspect of the main regularization schemes; what we 
can do is just to record them when they appear, and hope that in the future they 
can be better understood (Q). 


5 Some variations of the previous schemes 

The variations mentioned in the title are essentially of three kinds, described here¬ 
after in separate subsections. These variations will be mainly used in the applications 
of Parts II, III and IV; however, the hrst one will also be relevant for some subcases 
of the simple application proposed at the end of the present Part I (see subsections 
16.81 and 16.9p . 

^^One should probably look for their origin in some excessive idealization of the physical model 
(for example, one could try to describe in a more realistic manner the boundaries of the spatial 
domain; these are “hard” and “deterministic” in the present formulation, but could perhaps be 
replaced with “soft” or “stochastic” [43] boundaries). For a different type of boundary anomaly, 
related to vacuum effects for charged fermions, and for its cure taking into account back reaction 
effects, see [HI]. 
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5.1 The basic Hilbert space when 0 is an isolated point of the 

case of Neumann and periodic boundary conditions. In this subsection we 
are going to consider a variation of the framework developed in Sections [2] and |3] to 
deal with the cases where the fundamental operator A = —A + V acting on 
has its spectrum contained in [0,+cx3), with 0 an isolated point; in other terms, 
0 G cr{A) C {0} U +CX)), for some £ > 0 . In this case 0 is a proper eigenvalue, as 
it always occurs for isolated points of the spectrum. 

A standard approach employed in the physical literature to treat problems of this 
kind is to simply neglect the states of “zero energy”; see, e.g., [581160| IT0T| . According 
to the formulation considered in the present paper, this amounts to the following 
procedure: in place of we dehne the basic Hilbert space as the orthogonal 

complement in of the null eigenspace, that is 

Ll{Q) := (ker^)^ ( C L\Q) ) . (5.1) 

It should be noted that the restriction of A to is a selfadjoint, strictly positive 

operator in with spectrum contained in [e^,+cx)). In this situation, is 

the basic Hilbert space even from the viewpoint of held quantization (Q). 

Let us recall the dehnition (13.Ih . giving the integral kernel associated to a given 
operator on L^(H), and consider Eq.s (I3.17p . (I3.46p and (I3.53p for the Dirichlet, 
heat, cylinder and modihed cylinder kernels associated to A; if the latter operator 
is redehned as the restriction of —A + H to in the cited equations we should 

formally replace 5x, Sy with E^Sy where Eq is the orthogonal projection onto 

Lg(f2) (suitably extended to distributions, so that it can be applied to 5y^Ay)- With 
this modihcation, the expansions (I3.18p . (I3.47p . (I3.48p and (I3.53p for the kernels 
mentioned above hold again, using the eigenfunctions of A in Lg(f2) (Q). 

Typical conhgurations of the above type are those where A = —A, the spatial 
domain H is bounded and the held fulhlls either Neumann or periodic boundary 
conditions on dVt (@); indeed, in such cases the spectrum of A in L^(H) is purely 
discrete, 0 is an eigenvalue and ker^ is formed by the constant functions. Therefore 
Lq(H) , dehned via Eq. (15.ip . is formed by the functions with mean zero: 


Ll{n) 


{/ e L^a) 


dxfix) 



(5.2) 


this, we mean that the Fock space ^ of the quantized scalar field living in is the direct 
sum of all symmetrized tensor powers of Lg(n). 

^"‘As an example, in the case described by Eq. (IQ) the projection Eq onto Lq{Q) is given by 
Eof = f — voi(n) In /(^) (Vol(fi) is the volume of ft); the previous prescription makes sense 
as well for / = <5x and gives EqS^ = ■ 

^^As will be observed in subsection l5.3l the case of periodic boundaries would be more properly 
formulated in terms of a free field on a torus, but this is cause of no concern for the present 
considerations. 
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Let us mention that an analogous framework can be considered for slab configu¬ 
rations where x Neumann or periodic boundary conditions are 

prescribed on dfli x In these cases one works with the reduced operator Ai 
acting in L^(f2i); the latter must then be replaced with the Hilbert space 


:= (ker^i)-^ 


/ e L\n,) 



dxi /(xi) 



(5.3) 


and the basic Hilbert space for the full theory on H is ® . 

In the applications to be considered in the following, whenever 0 is an isolated point 
of the spectrum we will always assume that the fundamental operator A (resp. ^i) 
has been redehned so that it acts on the Hilbert space of Eq. fl5.ip (resp. 

L^(Hi) of Eq. (ESD). 

5.2 The case where 0 is in the continuous spectrum of .4.. Let us pass to the 
case where the fundamental operator A = —A+V is non-negative (cr^A) C [0, -l-cxo)), 
and 0 is in the continuous spectrum of A; then 0 has zero spectral measure and is 
a non isolated point of the spectrum (otherwise, it would be a proper eigenvalue). 
Here are two examples of this kind. To obtain them we consider the operator 
A := —A in L^(R'^), or the operator A := —A in L^(f2) where H is the half¬ 
space {x G R'^ I > 0}, and suitable boundary conditions, say Dirichlet, are 
specihed on dQ = = 0}. In these cases A has a complete orthonormal system 

of (improper) eigenfunctions (Ek)kgj^; with corresponding eigenvalues where: in 
the hrst case, K, = R'^ (with the Lebesgue measure dk), Juc(x) := (27r)“‘^/^e**^'^ , 
cuk := |k| ; in the second case, K. = (0,-I-cxd) x R*^"^ (again, with the Lebesgue 
meaure dk, J^k(x) := \/2(27r)“'^/^and, again, cUk := |k|. In 
both cases ct(. 4.) = [0, -|-cx)) and the spectrum is purely continuous. 

The case of A non-negative, with 0 in the continuous spectrum, cannot be treated 
with the approach of the previous subsection: there is no way to obtain a strictly 
positive operator by simply removing 0 from the spectrum. In a more physical 
language, infrared divergences cannot be simply ignored and we must devise a more 
sophisticated way to treat them, as we are currently doing for ultraviolet divergences. 
A natural approach to the problem is to represent .4, as a limit 


^:=“lim”A (5.4) 

where Ae is a selfadjoint operator, depending on a parameter £ G (0, £o) and such 
that the spectrum of Ae is contained in -|-cxd); the deformed operator Ae is used 
everywhere in place of A, and the limit e —)■ O’*" is performed only at the end, after 
zeta renormalization has been carried out. In particular, we dehne the deformed 
smeared field operator 

0- ;= {K-^Ae)-^/^$ (5.5) 
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and the corresponding deformed, regularized stress-energy tensor T^^(x) whose VEV 
is given by 

(0|f-(a;)|0) = 

~ f 2 ~ ~ 

/ y=x 

■ (O|0-(x)0-(2/)|O) ; (5.6) 

for the above VEV we have expression analogous to 03.31113.3,^ in terms of the 
deformed Dirichlet kernel 

y) := ^r"(x,y) = ^y) ) (5-7) 



with s = (n ± l)/2 . 

As mentioned above, in this generalized version of the local zeta regularization the 
limit e ^ O’*" must be considered only after the analytic continuation has been 
performed; in particular, we dehne 


{0\fy,y{x)\0)ren 


lim RP 

e^O+ 


/0|f™(a:)|0) . 


(5.8) 


The above renormalized VEV can be expressed in terms of the renormalized kernels 
and where 



lim RP 

£^0+ 



lim RP 

£—> 0 + 




, (5.9) 




X, y):= lim i?P 


p=0 


(x, y)) = lim RP 

e^0+ 


^a^^P>^(x,y)) ; 


these functions play a role very similar to the ones introduced in Eq. (I3.37p for 
a strictly positive A, and allow to express the renormalized VEV (15. 8 h as in Eq.s 

(1T39J1M1J). 

In the sequel we will write and T'^ (or T'^), respectively, for the heat and cylinder 
(or modihed cylinder) kernel associated to A^- Proceeding as in Section [3] we obtain, 
for 3?s sufficiently large. 


1 ^+oo 

7^!(x,y) = A:^(t;x,y) ; (5.10) 

r(s) Jo 

-I ^+oo 

D,^(x,y) = dtP-iT^(t;x,y) ; (5.11) 

the above formulas are the starting point to discuss the analytic continuation in s 
of the Dirichlet kernel D^, for any hxed e G (0,£o)- 
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Let us associate to A to the “undeformed” fundamental operator A the heat and 
cylinder kernels 

K(t;x, y) := y) , T(t; x, y) := e"*'^(x, y) , (5.12) 

as well as the modihed cylinder kernel 

f(t;x,y) := (x/J'V*^)(x,y) (5.13) 

(see subsection I3.10p : these can be represented as in Eq.s fl3.47p . fl3.48p and fl3.53p 
in terms of the eigenfunctions {Fk)k£ic and eigenvalues {oJk)k&K. of A. We stress that 
the functions 115.9p do not possess integral representations of the form 

fl5.1Up fl5.1ip with K^,T^ replaced by ib, T; in fact, the corresponding integrals for 
K,T are typically divergent. In the sequel, we will present a more subtle way to 
obtain from K or T (and T). 

Up to now we have not specihed any particular form for As- The following two 
choices seem to be natural: 

A,:=A + e\ (5.14) 

As := (VA + e^ . (5.15) 

The hrst one corresponds to the idea, widespread in the physical literature, to treat 
infrared divergences adding a small mass e [691 193] ; the second one is less familiar 
and is justihed by the considerations that follow. 

Assuming A^ to have either the form fl5.14p or fl5.15p . we readily obtain the following 
relations allowing to express the deformed kernels in terms of the analogous 

basic kernels K, T: 

Ae:=A + e^ ^ K%i]x,y) = K{t]x,y) ] (5.16) 

Ae-.= {\fA + ef ^ T^(t;x,y) = e"^*T(t;x,y) . (5.17) 

In particular, assuming the kernels K, T to be meromorphic functions of t, the above 
relations imply that the deformed kernels K^,T^ are meromorphic as well; in these 
cases, the deformed Dirichlet kernel admits integral representations analogous 
to fl3.94p fl3.95p . involving the Hankel contour i5- For example, one can write 

£>;(x,y) = ‘ jldte-‘F(t;x,y) . (5.18) 

Starting from the above representation we can derive explicit expressions for the 
renormalized functions more generally, for 

(«:2^+"Zi:(x,y)) (ue{-l,0,l,2,...}) , (5.19) 


-Di^l(x, y) := lim i?P 
2 £->0 
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that could be called “renormalized Dirichlet kernels” of order —n/2 . More precisely, 
let us assume the modified cylinder kernel T(t;x, y) associated to the fundamental 
operator M to be a meromorphic function of t in the neighborhood of the positive 
real half-axis, fulhlling the bound 

|f'(t; X, y)| ^ (7 for 3fJt ^-|-oo and some (7, a>0 . (5.20) 


Then, we obtain the following result, for n = —1, 0,1, 2,... (see Appendix IHj) : 


. (5.21) 


Let us remark that Eq. fl5.2ip has the same structure of Eq. fl3.98p . dealing with 
the Dirichlet kernel when A is strictly positive. In the stricly positive case, the 
cited result was derived rigorously (from Eq. fl3.97p L with no need to introduce a 
regulating parameter e; in the present framework, instead, it would be impossible 
to establish (I5.2ip without using the regulator e. 

One could derive results similar to Eq. fl5.2ip . involving the “renormalized deriva¬ 
tives”, e.g.. 



(5.22) 


where z, w are spatial variables; indeed, for n = —1, 0,1, 2,..., we have 


a,^Di1(x,y) = (-l)"+'r(u+2)Res(t-("+2)a,^f(t;x,y);0) (5.23) 


if dzwT{t ; X, y), as a function of t, fulhlls conditions of the form stipulated previously 
for T(l;x,y) (see, in particular, Eq. fl5.20p ). 


To conclude, let us discuss the pressure on the boundary in the present framework; 
the main point is the fact that, as in the case of strictly positive A, there are 
two possible prescriptions for the renormalized pressure. The hrst alternative is to 
introduce, at each point x G dfl, a deformed, regularized pressure with components 

PfW (0|f;y-(x)|0)rf(x) (!£{!,.... ci}) , (5.24) 

where n(x) = (n^(x)) is the outer unit normal to the boundary; we then dehne the 
renormalized pressure at x as 

pr”(x) := lim i?,P pf (x) . (5.25) 

e^0+ u=0 

The second alternative is to put 

p[^”(x) := f lini (0|f;j(x')|0)^e„^ n^(x) (5.26) 

yx'Gi2,x'^x / 

where (0|Tjj(x')|0)ren is dehned according to Eq. 05.81) at all interior points x' of the 
domain. The prescriptions 05.251) 05.26P can, in general, give different results: this 
happens, for example, in the case of a massless field on a wedge-shaped domain, to 
be discussed in Section 5 of Part II. 
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5.3 Some variations involving the spatial domain. In the literature, a scalar 
field fulfilling periodic boundary conditions is often considered. To give a rigorous 
description of this configuration, one should better give up to viewing as an open 
subset of and pass to a description in terms of tori. For example, it is customary 
to speak of a field on the hypercube (0, a)^ with periodic boundary conditions, where 
a > 0 is some given length. In the most precise description of this configuration, 
is not (0, a)^ but rather the d-dimensional torus Tf := R*^/ (aZ^ ~ (R/aZ)'^ (where 
Z is the set of integers, so that aZ = {..., —2a, —a, 0, a, 2a,...}) (Q). 

In some applications to be considered in the subsequent Parts II and III, the space 
domain is an open subset of R'^ but, in place of the Cartesian coordinates x = (x*), 
it is natural to use for it some curvilinear coordinates (g*)i=i,...,d = q; in these 
coordinates, the line element of R'^ will have the form 

d£^ = aij(q) dq^dq^ . (5-27) 

The above spatial coordinates induce a set of spacetime coordinates = q 

on R X where := t and the g®’s are as before; clearly, the spacetime line element 
ds^ = —dt^ + di'^ will have the form 

ds^ = 9^.u{q) dq^dq'' , 

5'oo := -1 , Qio = 9oi ■= 0 , 9ijiq) ■= ap(q) for i,j G {1, ...,d} . 

The analogue of Eq. fl2.16p in the coordinate system (g^) is 

f;, := (5-29) 

where is the covariant derivative induced by the (fiat) spacetime metric fl5.28p . 
In principle, the covariant derivative should appear in place of any derivative 
9^; however we are working with a scalar field and it is well-known that 

= dfj_f if / is a scalar function . (5.30) 

The situation is different when we consider second order derivatives, which explains 
the appearing of in fl5.29p . Let us recall that 

- T^^dxf (= Vj.^/) if / is a scalar function , (5.31) 


^®The considerations of subsection 15 .1 1 for the periodic case are easily rephrased in terms of the 
torus T^. The operator A := —A acting in L^(T^) has 0 as an eigenvalue, with kerXl formed by 
the constant functions; again, 0 is eliminated viewing A as an operator acting in 


LliTi) := {keTA)^ = {fGL^{Ti) 




dx/(x) = 0 
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where we are using the spacetime Christoffel symbols := + d^g^p — 

dpgpu)- The above computational rule is more efficiently implemented recalling that 
= t and using the space covariant derivatives Di corresponding to the line element 
fl5.27p : these rely on the Christoffel symbols := | a^^{diahj + djttih — dhCLij). From 
Eq. fl5.28p one easily infers that Fh = (for f, j, fc G {1, d}) are the only non¬ 
vanishing coefficients; so, Eq. fl5.3ip for a scalar function / on spacetime implies 

V.,/ = A,/ = a.,/-yj4/, 

Voi/ = dom-) = d,{dof) = V,o/ , Voo/ = doof . 


As a further variation of our schemes, we can stipulate the spatial domain to 
be an arbitrary d-dimensional Riemannian manifold, possibly non flat; in any co¬ 
ordinate system {q^)i=i^...^d = q of the Riemannian line element di'^ will have 
a representation of the form fl5.27p . (Of course the position O = Tf, considered 
at the beginning of this paragraph in relation to periodic boundary conditions, 
amounts to choosing for O a very simple, flat Riemannian manifold). Given any 
Riemennnian manifold O, we can associate to it the spacetime R x O equipped 
with the line element = —df^ -|- df'^; this takes the form fl5.28p in coordinates 

As a final variation, we can assume the space domain O to be an open subset of a 
Riemannian manifold and prescribe boundary condtions on dfl. 

Many results in Sections |2l |3]and|l]are readily adapted to the variations considered in 
this subsection for the space domain. An essential point in making these adaptations 
is to remember that, when an arbitrary coordinate system is employed, the second 
order derivatives of scalar functions must be intended in a covariant sense and the 
computational rules fl5.32p must be applied. 


6 The case of a massless field on the segment 

6.1 Introducing the problem for arbitrary boundary conditions. To con¬ 
clude the present Part I we present a simple application of our general formalism, 
namely a 1-dimensional model describing a massless scalar field living on a segment, 
with no background potential. This means that 

d = 1 , 12 = (0, a) (a > 0) , A= = 0) (6-1) 

where G (0,a) is the standard Cartesian coordinate (Q). The field is assumed 
to fulfill Dirichief, Neumann or periodic boundary conditions at the boundary 
dVt = {0}U{a}; we will deal with each one of these alternatives separately, in 
the subsequent subsections 16.6116.91 . 

^"^We could have used, in place of the rescaled spatial variable 

x\ := xVa G (0,1) , 

in terms of which, we would have obtained simpler expressions for the results reported in sequel. 
Yet, we choose not to employ this rescaled coordinate in order to make the comparison with 
kwnown results more straightforward. 
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In passing, we note that the setting described above is the d = 1 case both for the 
conhguration with two parallel hyperplanes and for the d-dimensional box (to be 
considered in Parts II and IV, respectively). 

Let us make some comparison with the previous literature about the Casimir effect 
for a scalar held on a segment. First of all, we wish to mention the book of Bordag 
et al. na (see Chapter 2) and the work by Fulling et al. [50]; these authors derive 
the total bulk energy, for several boundary conditions, using regularization meth¬ 
ods different from zeta approach. More precisely, [15] uses an exponential cut-off 
regularization followed by Abel-Plana resummation, while [50] employs essentially 
a point-splitting procedure. These authors also obtain the force acting on the end¬ 
points of the segment by differentiating the expression for the total energy with 
respect to the lenght of the segment (see the comments at the beginning of subsec¬ 
tion 14.4p . Let us also mention the paper [68] by Mamaev and Trunov deriving the 
stress-energy VEV for a massless scalar held on a segment in the case of periodic 
boundary conditions, via point-splitting regularization. 

In all cases where the present section has an intersection with [ISHSOlEH], our results 
are in agreement with these references (@). 


^®JVote added. After this manuscript was written, we became aware of a very recent paper 
by Mera and Fulling [TO] who consider a massive scalar field on a segment, regularized via an 
exponential cutoff, and compute the stress-energy VEV by the method of images, i.e., as a sum 
over infinitely many optical paths. (Even though different from zeta regularization, this approach 
is related to the cylinder kernel due to the cutoff structure.) The authors of [70] consider the zero 
cutoff limit, whose treatment can be understood as a renormalization, and also give a number of 
results on the zero mass limit. With a small addition to these results, presented hereafter, our 
renormalized stress-energy VEV can be shown to agree with the zero mass limit of the calculations 
in m- As an example let us consider the (0,0) component of the stress-energy VEV, i.e., the 
energy density, in the case of Neumann boundary conditions; in [70] this is expressed as the sum 
of a term diverging quadratically in the cutoff, and of other terms which have a finite limit when 
the cutoff is sent to zero. If one renormalizes removing the divergent term, one obtains from [70] 
that 

(0|Too(a;^)|0)r-en = 


m Ki{2mna) Ari(2m|a;^ -|- na|) 

27ra n tt la;^ -I- na| 

n—1 n— — oo ' ' 



^ +0O 

-) ^ Aro(2m|a;^-I-na|) 

n— — oc 


where m is the held mass and Kq,Ki are modihed Bessel functions of the second kind (see Eq.s 
(10) (15-17) and (51) of [70], here re-written with L = a, /3 = ^—1/4 and putting r = I = 0 to 
indicate the choice of Neumann boundary conditions). It is known that Ko(z) = — In z-I-0(1) and 
Ki{z) = 1/z -I- O(zlnz) for z —0+. So one expects that, for m —^ 0+, '>TT'J2n=i Ki(2mna)/n 
—l/n^)/(2a) = 7r^/(12a) and J2n=-oc + na\) —^ 0; these two facts are 

in fact established in m- Let us add to these results the remark that, for m —>■ O’*", 
+na\)/\x'^+na\ -)■ (1/2) l/(^^ + = 7’‘^/(2a^sm^(^a:^)) (for 

the computation of the last series by contour integral methods see [59], p. 268, Eq. 4.9-4). In view 
of these facts 


lim (0|roo(a:i)|0) ren — 
771— )-0 + 




24o2 ^ 2a2sinV2la:i) ’ 

^ a ’ 


in agreement with Eq. (16.401) of the present work. Let us emphasize that, differently from [70] . 
our zeta approach gives the renormalized stress-energy VEV by mere analytic continuation, with 
no need to rimove divergent terms. 
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6.2 Cylinder and Dirichlet kernels. For any one of the previously mentioned 
boundary conditions, we perform our analysis in the manner explained hereafter. 
First of all, we determine explicitly the cylinder kernel T(t;a;^,?/^) associated to 
the fundamental operator A; to this purpose we consider a complete orthonormal 
set of eigenfunctions {Fk)keK for ^ with eigenvalues {ul)k£K- The label set /C is 
countable and Jj^dk means eigenfunction expansion fl3.48p for the 

cylinder kernel reads 


T{i-x\y^) = • (6.2) 

fce/c 


Once the cylinder kernel has been computed explicitly by evaluating the above sum, 
we can proceed to determine the modihed cylinder kernel T(t; y^) as the primitive 

of — T(t; vanishing for t — )■ +cxd (see Eq. (I3.55p ). 

In the next subsections T and T will be computed explicitly, for several kinds of 
boundary conditions. In all cases, the cylinder kernel T and the spatial derivatives 
of the modihed cylinder kernel T will be found to have meromorphic extensions in 
the t variable to an open complex neighborhood of [0, +C)o), with possible poles only 
at t = 0, vanishing exponentially for +cx3; thus, the framework of subsection 

13.151 can be applied straightforwardly. 

To be more precise: for y^ , T,T and their derivatives have analytic extensions 
in t to a neighborhood of [0,+oo). When evaluated on the diagonal y^ = x^, the 
modihed cylinder kernel T is found to have a logarithmic singularity in t = 0, while 
the cylinder kernel T and the spatial derivatives of both T and T are meromorphic 
in a neighborhood of [0, +cxo) with only a pole singularity in t = 0 . Because of this, 
one can resort to Eq.s fl3.95p and fl3.97p to obtain the analytic continuation of the 
Dirichlet kernel and of its derivatives, required in order to determine the regularized 
VEV of the stress-energy tensor; explicitly, we have 


Du=i{x\y^) 


3 —27r(u—1) 


F(2- 


-u 


27ii 




(6.3) 


dzwDu+x {x , y ') 


3 —i7r(u+l) 


F(l- 


-u 


2Tii 




(6.4) 


for z,w E {x^,y^} . 


In order to obtain the analytic continuations at m = 0 of the above functions, one 
can simply set m = 0 in the expressions on the right-hand sides of Eq.s fl6.3p fl6.4p 
and explicitly evaluate the remaining integrals along the Hankel contour via the 
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residue theorem (j^; as indicated in subsection 13.151 this gives 


D_i{x\y^) =-Res{i 


(6.5) 


d^y,Di{x^,y^) = Res t ;0 ioi z,w e{x^ ,y^} (6.6) 

2 yl—^1 \ yl—^1 / 

(use Eq. fl3.96p and the analogue of Eq. fl3.98p for the derivatives of Dg). 


Before moving on, let us mention that analogous considerations can be made con¬ 
cerning the traces Tr^“*, T(t) (see Eq. fl3.22p and Eq.s fl3.56p fl3.57p . respectively). 
Indeed, we can compute the cylinder trace T(t) according to Eq. fl3.58p that in the 
present case reads 

T(t) = (6.7) 

keic 

By explicit evaluation of the above sum for the boundary conditions considered in 
the sequel, it becomes apparent that T(t) possesses the same features of its local 
counterpart. Thus, we can resort again to the general framework of subsection 13.151 
to obtain the analytic continuation of Tr^“^; in particular, due to Eq. fl3.99p . the 
continuation at s = —1/2 is 


Tiv 11/2 = -Res(r'^r(t);oi , 


( 6 . 8 ) 


6.3 The stress-energy tensor. We can now determine explicitly the renormal¬ 
ized VEV of the stress-energy tensor; in fact, since no singularity arises, Eq.s (13.391 
I3.4ip and fl3.42p imply 


(0|foo(x^)|0)^en = + (^-^d^iyiDi{x^,y^) 


(6.9) 


(0|foi(x^)|0)^en = (0|fio(x^)|0)ren = 0 , 
{0|fii(x^)|0)ren = 

^'^D_i{x^,y^) + ^ d^iyiDi{x^,y^) - ^ d^i^iDi{x^,y^) 


( 6 . 10 ) 

( 6 . 11 ) 


In the following, the scheme outlined above will be illustrated in detail, as an exam¬ 
ple, for the case of Dirichlet boundary conditions. For the other boundary conditions 


^®Notice that the analogoue of Eq. (16.31) in terms of T is not so simple and straightforward to 
employ, due to the logarithmic behaviour of the modified cylinder kernel near i = 0; on the other 
hand, the analogue of (|6.4I) in terms of T has a singularity in the gamma function for u = 0. Thus, 
for the computations in which we are interested, there is no better strategy than using Eq. (EH) 
with T and Eq. (16.411 with T ; this also explains why, in the sequel, we will frequently refer to both 
kernels. 
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we will be more synthetic but, in any case, we will always report the expressions 
for r(t;x,y), r(t;x,y) and (0|T^i,|0)ren; in particular recall the considerations of 
subsection 12.51 and note that in this case Eq. fl2.24l) gives 


6 = 0 . 


( 6 . 12 ) 


6.4 The total energy. Since no singularity appears, we can use the general pre¬ 
scription (14.1111 ) the latter, along with Eq. (16.Sp . allow us to derive an explicit 
expression for the bulk energy, for any one of the several boundary conditions to be 
considered in the following. More precisely, we have 

i Res(^r2T(t) ;0) (6.13) 

where T(t) is the cylinder trace of Eq. (I3.57p . 

In passing, let us remark that the renormalized boundary energy always van¬ 
ishes identically in the cases considered, due to the prescribed boundary conditions 
(indeed, the same statement can be made for the regularized version see Eq. 

dm)). 

We also mention the following fact: by direct comparison of the results reported in 
subsections ESES] it appears that the results derived using Eq. (I6.13P could as well 
be deduced integrating over (0,a) the conformal part of the renormalized energy 
density (0|Too|0)r.en • On the contrary, the non-conformal part of the latter appears 
to diverge in a non-integrable manner near the end-points x = 0 and x = a. 


6.5 The boundary forces. Let us remark that, since the boundary is zero¬ 
dimensional, the nominal “pressure” on the boundary points = 0, = a does in 

fact coincide with the force on these points; because of this, we adopt the notation 

Fren{x^) = p”'^”(x^) for x^ = 0, o . (6.14) 


For all the (non periodic) boundary conditions to be analysed in the following sub¬ 
sections, there are in principle two dehnitions of the renormalized boundary forces; 
these descend from the two alternatives pointed out in the general discussion on 
pressure of subsection 14.21 

Let us indicate with n^{x^) the unit “outer normal” at the points on the boundary, 
so that ?7,^(0) = —1 and n^(a) = 1. The hrst dehnition reads 


(x^) := (0|f“(x^)|0) 




4i=0 


(6.15) 


(see Eq. (I4.22p and notice that the prescription of taking the regular part is superflu¬ 
ous, since no sigularity arises; namely, we hrst compute the regularized stress-energy 
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tensor at the boundary point and then we analytically continue at u = 0). The 
second alternative is to define (see Eq. fl4.23p l 


X : = 


lim 

x'^G{0,a), x'^- 


(0|Tii(x'^)|0),.en ) n^(x 


(6.16) 


(i.e., we first renormalize at inner points of the interval (0, a), and then move towards 
the boundary). As a matter of fact, for all boundary conditions considered in the 
next subsections, the equivalence between fl6.15l) and fl6.16p will be checked by direct 
computation. 

6.6 Dirichlet boundary conditions. As a first example, let us consider the case 
where the field fulfills Dirichlet conditions at both the end points of the segment 
(0, a), that is 

x^) = 0 for f G R, and x^ = 0 or x^ = a . (6.17) 


A complete orthonormal set of eigenfunctions {Fk)keic for A and the related eigen¬ 
values {ul)keK are 


.(x^) := \ - sin(/cx^), ul := k‘^ 

V a 


for /c G /C = 


nvr 


n = 1,2,3,... \. (6.18) 


The expansion 06.21) for the cylinder kernel associated to A reads 


+ CXD 


rp/, 1 In ^ -iHLt . fnTT , /nvr 
I [i] X ,y ) = - y e “ sm — x sm — y 
a \ a / V a 

n=l 


(6.19) 


re-writing the trigonometric functions in terms of complex exponentials, the right- 
hand side of the above equation reduces to a sum of four geometric series, which 
can be explicitly evaluated. The final result is 


T{i;x\y^) 


1 cos(^(x^—j/^)) — e A 

2a cosh(^t)—cos(^(xi—j/i)) 


cos(^(x^-|- 2 /^)) — e A 
cosh(ft)-cos(^(xi + 2/i)) 


( 6 . 20 ) 


the same expression is also reported, e.g., in [13 |50], but therein it is not used to 
compute the full, renormalized stress-energy VEV. 

Expressing the hyperbolic functions in terms of exponentials, we easily obtain the 
primitive of T which vanishes exponentially for t —?• -fcxo, that is —T; in conclusion 


T{i-,x\y^) = 


27r 


TT 


In 1 — 2e “‘cos ( —{x^—y )) + e “ ‘ M- 


Infl —2e “‘cos (x‘-|-|/‘)j-I-e “ 


( 6 . 21 ) 
















Both T and the space derivatives of T have meromorphic extensions in t to a complex 
neighborhood of [0, +cxo), with poles only at t = 0; so, we can employ Eq.s flh.Sp fl6.6p 
to obtain from them the renormalized Dirichlet kernel and its spatial derivatives. 
For example, since 


1 

vrt 


12a2 sin^(- 

^ a ' 


i + 


ni-x\y^) 

7r3(15(2 + cos(^a:i))-sin^(f a;!)) 


720a4 sin^(^a;i) 


F + 0(F) , 


( 6 . 22 ) 


for t —)■ 0, evalnating explicitly the residne in Eq. (I6.5p . it follows 



TT 3-sin2(^a;i) 
120^ sin^(-a;^) 

k a ' 


(6.23) 


Proceeding similarly for the derivatives of the Dirichlet kernel, and then nsing Eq.s 
(16.9116.11]) . one obtains the following expression for the renormalized VEV of the 
stress-energy tensor: 


(0|f^^(a;^)|0)ren 

A-=^ 

■ 24a2 ’ 


= A 
B(x^) : = 


-1 0 
0 -1 

TT 1 


+ ^B(x^) 


1 0 


sm 


2fn 


X^) 


0 0 
for x^ e (0, a) . 


(6.24) 


Let ns now discnss the renormalized bnlk energy. To this pnrpose, we first note that 
the expansion (I3.58p for the cylinder trace gives 


+ 0O 


T(t) = Z‘^ 

n=l 


_ riTT ^ 




1 ’ 


then, using prescription f l6.13p . we readily infer 

TT 

^ren _ _ 


24a 


(6.25) 


(6.26) 


In conclnsion, let ns consider the bonndary forces; it is easily seen that both defini¬ 
tions (16.151) and (16.161) give (with A as in Eq. (16.241) ) 


Pren(O) — A , Fren{,o) — —A . 


(6.27) 


6.7 Dirichlet-Neumann boundary conditions. Let ns now consider the case 
where Dirichlet and Nenmann bonndary conditions are respectively prescribed at 
the two end points of the segment (0, a): we assnme 

0(t, 0) = 0 , dxi4>{t, a) = 0 for t G R . (6.28) 
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In this case, a complete orthonormal set of eigenfnnctions (-Ffc)fceK; for ^ and the 
related eigenvalnes {ojDkeic are described by 


Fi 


{x^)-.= sin{kx^), ul ■.= for /cG/C = |^n+- j — n = 0,1, 2, ...j-. (6.29) 


Using the expansion fl6.2p . we can determine the cylinder kernel T and then obtain 
the modified kernel T as minus the primitive of T, vanishing for t —>■ +cxd; the hnal 
results are 


T{i-,x\y^)=- 

a 


sinh(^ t) cos(^(xi- 2 /i)) sinh(^ t) cos(^(a:^+y)) 


cosh(^ t) — cos(^(a:^— 1 /^)) cosh(^ t) — cos(^(a:^+?/^)) 


, (6.30) 


Tit;x\y^) = 


1 


pYcos(^(a:^-|/^)) + cosh(^t) 


In 


/ cos(^(3;H|/^)) + cosh(^t) Y 
Vcos(^(a;i+2/i))-cosh(^t)y_ ’ 


(6.31) 


Vcos(^(a:i-|/i))-cosh(|^t) 

Using the above expressions along with Eq.s (16.51) . (16.6p and (16.9116.lip , one obtains 
the renormalized VEV of the stress-energy tensor: 


(0|f^^(x^)|0), 


= A 


A : = 


TT 


48a2 ’ 


/j.,iy=0,l 


1 0 
0 1 


+ es(x') 


1 0 
0 0 


TT cos(-a:^ 


(6.32) 


TTT-27y—iT fo’^ ^ (O’ 

2a2 Sin^ ^ ^ 

\ n / 


Next, we derive the cylinder trace using again the expansion fl3.58p : 

+0O 

n=0 


+ 0O TT , 


eo. 


-‘-1 


(6.33) 


Now prescription fl6.13p allows us to obtain the renormalized total bulk energy: 

71 


^ren _ 


48a 


(6.34) 


Concerning the boundary forces, also in this case dehnitions fl6.15p fl6.16p agree and 
give 

Fren{0) = , -^ren(a) = A (6.35) 

where A is as in Eq. fl6.32p : notice, in particular, that the above expressions have 
the opposite sign with respect to the ones of Eq. fl6.27|) . corresponding the case of 
Dirichlet boundary conditions. 
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6.8 Neumann boundary conditions. We are now going to study the case where 


= 0 for t G R, and = 0 or = tTq . (6.36) 

In this case, according to the considerations of subsection 15.11 the Hilbert space 
L^(0, a) has to be replaced with the space hg(0, a) of square integrable functions on 
(0,a) with mean zero (see Eq. fl5.2p ): in this space, a complete orthonormal set of 
eigenfunctions {Fk)k£K for ^ = —dx^x^ and the corresponding eigenvalues {Uk)k&K. 
are given by 


Efc(x^) 



cos(fcx^), 



for k & K, = 




(6.37) 


The cylinder kernel associated to A can be evaluated according to Eq. fl6.2p to 
obtain 


T{t;x\y^) 


1 r cos(^(x^ — 1 /^)) — e ‘ cos(^(x^+|/^)) — e ‘ 

2a cosht — cos(^(x^ — ?/^)) cosht — cos(^(xi+j/i)) 


(6.38) 


while for the modihed cylinder kernel, computed as minus the primitive of T, we 
obtain 




1 


In 1 


_ TT 1 

2e “ cos 



+ 


+ In 


1 - 2e 


—-t 
“ COS 




(6.39) 


Resorting once more to Eq.s fl6.5p . fl6.6p and fl6.91l6.lip . the renormalized VEV of 
the stress-energy tensor is found to be 


(0|T),^(x^)|0), 


= A 




-1 

0 


^B(x^: 


(6.40) 


where A and B(x^) are dehned as in Eq. fl6.24p . 

In the case under analysis the spectrum of A coincides with the one obtained for 
Dirichlet boundary conditions (compare Eq.s flb.lSp fl6.37p h therefore, the cylinder 
trace is again given by Eq. fl6.25p , and we derive the same renormalized bulk energy 
as in Eq. (I6.26p : 

vr 

£jren _ _ __ 

24a ■ 

Finally, both dehnitions fl6.15p and fl6.16p give for the boundary forces the same 
results as in the case of Dirichlet boundary conditions (see Eq. fl6.27p ): 


-^ren(O) — A , Fren{(^) — —A 


(6.41) 


(again, A is as in Eq. fl6.24p h 
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6.9 Periodic boundary conditions. The last case we consider for the segment 
conhguration is the one where the held satishes periodic boundary conditions: 

4>{t,0) = (p{t,a) , dxi(p{t,0) = dxi(p{t,a) for t G R . (6.42) 

As explained in subsection l5.31 this case would be more properly formulated in terms 
of a free scalar held on the 1-dimensional torus := R/(aZ) . Besides, similarly 
to the case of Neumann boundary conditions, recall that the basic Hilbert space is 
Lq(T^) = {/gL^(T^) I dx^f{x^) = 0} (see subsection 15.11 and the footnote 1261 of 
page[S2])- In this space a complete orthonormal set of eigenfunctions {Fk)k£K for 
with the corresponding eigenvalues is 


Fu(x^) : = 


Ak 


2 / 2 
Uk ■= k 


for keJC = < ± 


2n7r 


n = 1,2,3,...} . (6.43) 


Let us pass to determine the cylinder and modihed cylinder kernel associated to A; 
using the same methods of the previous subsections, we obtain 


r(t;x ,y ) = 


cos{^{x^ — y^)) — e 
a [cosh(^t) - cos(^(a;i-?/i))] 


27r 


_27r I / 27r 


T{i-,x^,y^) = -In 1 — 2e “ mos ( — {x^ — y^)] + e 


V a 


_ TT ^ 


(6.44) 


(6.45) 


(the same expression for T is also reported, e.g., in |17], again for other purposes). 
Eq.s fl6.5p . fl6.6p and (16.9116.lip , yield the following expression for the renormalized 
VEV of the stress-energy tensor: 

'-1 0 
0 -1 


(0|f^,.(a:^)|0), 




TT 


(6.46) 


Let us stress that the above results respects the invariance under translations x^ i—)■ 
x^+a (for any a G R) of the given conhguration, since it does not depend explicitly 
on the spatial coordinate x^. 

To conclude, we discuss the renormalized bulk energy. We hrst note that expansion 
fl3.58p for the cylinder trace yields, in the present case. 


-1 


+ CX> 


H-cxd 


r(‘) = E + E 


‘ = 2E- 


^TITV ^ 


\n=—oo n=l 


n=l 


e a ^ — 1 


then, using once more prescription fl6.13p . we obtain for the bulk energy 

TT 

^ren _ _ __ 

Qa 


(6.47) 


(6.48) 
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A Appendix. On the form (l2.11lj for the stress- 
energy tensor 


Eq. (12.lip is the quantized version of a classical formula for the stress-energy tensor 
[121 [13 UHl ESli which we review here for completeness. Following Section [3 we 
refer to Minkowski spacetime; after an inertial frame has been chosen, the latter 
is identified with = R x R'^ 3 x = (x^) = (t,x). We confine the attention 

to a subset of the form R x hi, where C R'^ is a spatial domain; let us consider 
a classical scalar field 0 on R x described by an arbitrary Lagrangian density 
C = C{(f),d(j),x). The associated canonical stress-energy tensor is 

djC 

and fulhlls 

^ ^A.2) 

along the solutions of the held equations. If S is any spacelike hypersurface with 
normal unit vector and volume element dv, we dehne the canonical momentum 

F“^(E) := j . (A.3) 

For simplicity we rescrict the attention to the case 

S = {t} X , (A.4) 

for a hxed f G R, writing for the corresponding canonical momentum. In 

this case we can take {N^) = (1, 0,..., 0) and dv corresponds to integration on 
n with respect to the usual volume element dx, so 

pr{t)-.= [ d^TsrM- (a.5) 

Jn 

Writing x for (f,x), we have dP^°''^ / dt{t) = J^dxdoTQ^^{x) = dx (—0°Tg“")(x) = 
— d^T^l^){x) = dx [0*T“"'(x) -|- (0i.£)(0(x), cI0(x), x)], i.e., by the 
d-dimensional divergence theorem, 

—j— (t) = / dx (0,^£)(0(x), 00(x), x)-f / daipi) rd{-x) {x) . (A.6) 

dt Jn Jqq 

Here (and in the sequel) n(x) = (n*(x)) is the outer unit vector in R*^ normal to 
the boundary dQ at x, and da is the (d —l)-dimensional area element (Q). For a 

^°Obviously enough, if il is unbounded we intend fg^da(x) njx) := fgQ^da£(x) n}(x) 

where ili C 112 C tla C ... are bounded domains such that = H. 
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number of reasons briefly reviewed in the sequel, it is customary (see, e.g., mm) 
to consider an “improved stress-energy tensor” of the form 


rji _ rpcan I oA 

fiiy • ^ Xfiu 1 

where F\^y is a covariant tensor of rank 3 such that 

Fxfiu F^xu ■ 

Condition (lA.Sp implies d^{d^Fx^u) = 0, thus ensuring 




(A.7) 


(A.8) 


(A.9) 


We can give dehnitions similar to flA.3p and flA.Sp using the improved stress-energy 
tensor; in particular, for each f G R, we dehne the “improved momentum” 


P^{t) := / dxTo^(t,x) 
Jn 


(A.IO) 


We claim that 


P^{t) = + / da(x) n*(x) FioX^^^) , 

Jan 


(A.ll) 

where da and '^•(x) have the same meaning as before. To prove this, note that 


-Ou 


rpcan 

-^Oiy 


= d^Fxo^ = d^Foo, + = d^F, 


iOu 


(A.12) 


(^oor “ 0 due to Eq. IIA-8[) E thus P^t) — + /,jdxt?’F,or(t,x), and the 

d-dimensional divergence theorem yields Eq. flA.llll . 

In many cases of interest, the boundary term in Eq. flA.llll is zero. In particular, 
this happens if hi = R'^ and FiQi,{t,x.) vanishes rapidly for x —)■ oo. In the case of 
a bounded domain, the boundary term can be zero if depends suitably on the 
held (j) and the latter fulhlls appropriate conditions on dQ. 

Whether or not the boundary term in Eq. (1A.11I1 vanishes, using Eq. (1A.9I1 we prove 
that the improved momentum evolves according to the analogue of Eq. (lA.bp . i.e.. 


= [ dxc {di,C){(j){x),d(j){x),x) + [ da{x.) n\x.) Ti^{x) 
Jn Jan 

The improved stress-energy tensor is symmetric if and only if 


d\Fx^u - Fa.;.) = - (T; 


can 

fiiy 


rpcan\ 

Ufl ) 


(A. 13) 


(A. 14) 


When T™"" is not symmetric and it can be found a rank 3 tensor Fxfj,^ fulhlling 
conditions flA.Sp and flA.ldp . the symmetry of the improved stress-energy tensor 
(IA.7P is itself a good reason to consider this object. 
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There are reasons to consider the improved stress-energy tensor even in the case 
when is itself symmetric (of conrse, in this case Eq. flA.ldj) reqnires that 
be symmetric in /i and v). One of these reasons has been pointed ont by 
Callan et ah [19]; in few words, after qnantization the divergences of the improved 
tensor can happen to be softer than the divergences of the canonical one, especially 
in pertnrbative renormalization. 

In this paper we are interested in a held theory governed by the eqnation 0 = 


{—dtt -l- A — l/)0 = {d^d^ — V)(j) which arises from the Lagrangian 

. (A.15) 

The corresponding canonical stress-energy tensor is 

I;T = 1 + V4>^) ; (A.16) 

this is symmetric and fnlhlls (along solntions of the held eqnations) 

= (A.17) 

Condition (1A.8I) is satished by the tensor 

Fxi^u ■= , (A.18) 

where ^ is a real parameter. In this case 

d^Fx^,y = - r]f,^d^dx)(j)'^ ; (A.19) 

this tensor is symmetric in fi and u, so is symmetric as well. The derivatives of 
0^ in Eq. flA.19p can be re-expressed nsing the held eqnation = V(j), yielding 

d^Fx^. = -2e d,(j) + 2e dx(j) + Vcj)^) - 2^ 0 . (A.20) 

Tims the improved stress-energy tensor (lA.7p takes the form 

= (1-2O0;.00.0 - Q - 2^^ r]^,{d^ct>dxct> + Vcfa^) - , (A.21) 

of which (I2.11|l is a natural quantization. 


Let us recall that the momentum corresponding to the improved tensor is 
related to the canonical one via Eq. (lA.lip : in the present framework where Fx^u 
is given by Eq. flA.18p . the boundary term in Eq. flA.lip vanishes under Dirichlet 
boundary conditions (0(f,x) = 0 for all x G dQ)] this term vanishes as well if 
n = R'^ and 0(f, x), 0a0(C x) vanish rapidly for x —)■ cx). 


65 



















In the case V = 0, the above improved tensor flA.2ip allows another interpretation: 
this is the functional derivative with respect to the metric of the action describing 
a scalar field that interacts with a gravitational field via the scalar curvature, in the 
limit where 0 is small and the metric is the Minkowski metric 77 ^ 1 , plus a perturbation 
of the second order in <p. Concerning this, see naiisiEi] and |3Sl EZ] • 

Again for V = 0, the action of the field coupled to gravity is conformally invariant 
for ^ = (d- l)/(4d) [TnT] . 

B Appendix. Smoothness properties of some in¬ 
tegral kernels 

We refer to an operator A in with features described by Eq.s fl3.7p . fl3.8p or 

fl3.9p (recall that the situation fl3.7p is more general than fl3.8p . and fl3.8p is more 
general than fl3.9p . In the present appendix we present a number of results which are 
useful in relation to several integral kernels associated to A . Proving these results 
would require a heavy use of functional analysis, which is not among the purposes 
of the present work; therefore, hereafter we just sketch the basic ideas that will be 
presented with more details elsewhere [4T] . 

General results for case (13.711 . As in the cited equation we assume A to be 
any strictly positive, selfadjoint operator in A basic step towards our goal 

consists in associating a scale of Hilbert spaces 'H'’ to the powers A^ for r G R (see 
[H] for a similar construction); speaking somehow loosely, we can describe 'H'" as the 
Hilbert space of generalized functions / : H —)■ C such that A''f G -h^(hl), equipped 
with the inner product 

(/Is). := (X-'VlX-'^g) (B.l) 

(as usual, ( | ) denotes the inner product of L^(r 2 )) which induces the norm 

ll/IG= \/W7 = IM’-'Vll (B.2) 

Let be any complete orthonormal set of (generalized) eigenfunctions 

^^Let us sketch the precise definition of W, given in [41]. For any real r we consider in L^(f2) 
the dense linear subspace O'" on which is well defined according to the general framework 
for functional calculus of selfadjoint operators in Hilbert spaces (see, e.g., [HI]); if {Pk)kinK is a 
(generalized) complete orthonornmal set of eigenfunctions of A with eigenvalues (wfc > £ > 0), 
then is formed by the functions / G L^(H) such that Ji^dkojl'^\{Fk\f)\^ < + 00 . We introduce 
on D'' an inner product ( | )r following Eq. and define 'H’" to be the completion of I?’' with 

respect to ( | jy ; for r 2 ^ ri one finds C a fact implying C Of course 

coincides with the identity operator 1 of L^(f2); thus = lA' = and ( | )o is the usual 

inner product ( | ). For any r ^ 0, is itself complete, so TL^ = C for r < 0 the 

completion TL^ is larger than D^, and even of L'^{Vl) = Ti,^. If r is sufficiently large, T-L~^ contains 
elements which cannot even be interpreted as ordinary functions on H; for example, as shown later 
in this appendix, the Dirac delta at any point x G D can be interpreted as an element of 
for all r > d/2 . 
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of A with eigenvalues {ul)k£K ^ > 0 for all /c G /C); then W is made of the 

generalized function f on Q such that 


/ dkut\{F,\f)\^<+oo; 
Jk 

moreover, for any pair of functions f,g E W, there holds 


(B.3) 


(/l9)r = [ dk uf (Ft\S) (fils) . (B.4) 

Jk 

Of cours^'H° is the usual space -h^(O), and it can be proved that 1-U"^ for 

f '2 ^ "Ti Moreover, from the characterization flB.Sp of l-F and from the identity 
A^Fk = ojl^Fk it follows that 

A~^ maps continuously into for r, r' G R such that r' — r < 23fJs ; (B.5) 


similarly, one proves that 

e~^'^ map continuously into for all r, r' G R . (B. 6 ) 

Next, consider the inner product of for all r G R this can be extended to a 

continuous, sesquilinear map 


n } :n-^ xn^ , (B.7) 

which allows, in particular, to infer the isomorphic identihcation (?{^)' = 'H~^. 

Results for case (13.8p . In accordance with the cited equation, we now assume 
A = —A + V on an open subset hi of R'^ with given boundary conditions, 1/ : hi —)■ R 
a C°° potential; all these ingredients are chosen so as to ensure the strict positivity 
of A. In this case it can be proved that 

^ C\Q,) for j G N, r G R with r > ^ + j , (B.8) 

where C^{Q) carries the topology of uniform convergence of the derivatives up to 
order j on compact subsets of hi; this is induced by the family of seminorms 

PjkU) {K C kl compact) (B.9) 

|o|s£j, yc£K 

Given any pair of topological vector spaces X,y, we say that X is continuously embedded in 
y, and we write X ^ y, ii X is a linear subspace of y and the identity map from X into y is 
continuous. So, for example, ^ ^ F{n) ^ ^ TF^ if ^4 ^ ra ^ 0 ^ r 2 ^ ri. 
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(B^IO) 


(0). By duality, Eq. flB.Sp implies (under the same assumptions on j and r) 

{c^{n))' ^ {vry = n-'". 

Now, let j G N and x G hi; the prescription 


(<5x,/) := /(x) 


(B.ll) 


makes sense for all / G (hi) and allows to interpret the Dirac delta function as 
an element of the dual space (C'-^ (r2))'. Moreover, it can be shown that the mapping 
h : X G D t—)■ ^x; G (C-^ (D))' is itself of class (0); of course, this fact and Eq.s 
fIB.lOD imply that 


5 : D — )■ , X !-)■ ( 5 x is CMf r G R, j G N and r > j . 

Finally, let j G N and suppose that 

B : 'H-('^/2+i2+»?) '^d/2+n+r] jg lingaj; and continuous 

for some rj > 0 and all ji , j 2 G N with ji + j 2 ^ j ; 

we claim that, in this case. 


(B.12) 


(B.13) 


the kernel D x D —)■ C, (x, y) i—)■ B(x,y) := (5x|i3(5y) 
is (well dehned and) of class C^(Q x Q). 


(B.14) 


In fact, each derivative dB( , ) of order ^ j involves (in an arbitrary order) 
operations of derivation with respect to x* and jSi operations of derivation with 
respect to y* (z = 1,..., d) where ji := cti + ... + ad and j 2 := /3i + ... + fdd are such 

^^Let us sketch the derivation of Eq. (IB. 81) . Due to well-known results on elliptic operators [25] . 
one has 

T-T ^ H[^yQ) for all r 0 , 

where is the standard local Sobolev space of order r (that is, the space of functions 

/ : D —>■ C such that (1 — A)’’/^(y/) € L^(fi) for all smooth, compactly supported functions 
y : D —>■ C). On the other hand, the usual Sobolev imbedding theorems give 

Hloci^) ^ C^{U) for all r S R, j S N with r > j . 

Summing up, the embeddings discussed in this footnote yield the thesis (IB. 81) . 

^^For example, setting j = 1 , one finds that the map x i—>■ Jx G (C'^(B))' is with derivatives 
{diS)x such that 

((d*5)x,/) = dj(x) (f G {1, ...,d}) 

for all / G C'i(D). 
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that ji + j 2 ^ j] any such derivative exists and is continuous on x with the 
explicit expression 5i3(x, y) = {{d°‘S)y:\B {d^6)y) 

Applications to the Dirichiet, heat and cylinder kernels: regularity results. Let 
s G C, j G N; if3fJs > d/2 + j/2, using Eq. fIB.Sp one easily infers that the 
operator B := fulhlls the condition (IB.lSp (with rj = lR.s — d/2 — j/2); thus 
flB.ldp holds. In conclusion, we have the following result for the Dirichiet kernel 
^.(x,y) := 

Da G {Vt X D) for s G C, j G N and > - + ^ . (B.15) 

Similarly, for any t > 0 and any j G N, due to flB.dp the operators B = or 
B = fulhll the condition fIB.lSp . implying Eq fIB.ldp . So, the heat and cylinder 

kernels iL(t;x,y) := T(t;x,y) := ((5x|e“‘^5y) are of class 0 in x,y, 

for all j G N; in conclusion, these kernels are in x, y; 

Ls:(t; , ),T(t; , ) G X D) for each t > 0 . (B.16) 

Results for case (13.911 . We now consider the case described by the cited equation. 
Thus A = —A + V on bounded domain D with C°° boundary dQ, on which Dirichiet 
conditions are imposed; the potential V is in C'°°(D) and E(x) ^ 0 for all x G D 
(recall that, due to these assumptions, A is selfadjoint and strictly positive). 
Smoothness up to the boundary. In the setting described above, one can strengthen 
the results discussed in the previous paragraph of this appendix replacing system¬ 
atically the space (7'^ (D) with C-^/D); this space consists of the functions which are 
continuous with their derivatives up to order j on the closure D = D U dfl, and it is 
equipped with the norm 


ll/llc. := max_|a7(x)|. (B.17) 

|a|sSi,xeO 


Similarly, the space C^{fl x D) can be replaced with C^{Q x D), endowed with the 
norm 

WfWci ■= max _|5x^^/(x,y)l • (B.18) 

|a| + l/3Ki> x,yeO 

In the present case Eq. fIB.Sp has the stronger version 

^ for j G N, r G R with r > ^ + j ; (B.19) 

^^note that, due to (IB.121) . the map y i—>■ {d^S)y is continuous from Q to due 

to (|B.13ll . y !->• B{d^6)y is continuous from tl to finally, due to (IB.121) the map x >->• 

(d“(5)x is continuous from ft to ^ and the sesquilinear form ( | ) is continuous on 

'^-(d/2+jl+i)) X '}{d,/2+ji+ri ^ 
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similarly, Eq. flB.14p (with the assumptions flB.lSp on B) holds with C^{VL x hi) 
replaced by x hi) (Q); due to these facts, Eq.s flB.lsp and (IB.16P for the 

Dirichlet, heat and cylinder kernels hold in the stronger versions 


Ds e X fi) for s G C, j G N with sRs > ^ + | ; (B.20) 

K(t; , ),T(t; , ) G X H) for each t > 0 . (B.21) 

Estimates for the eigenfunctions and eigenvalues of A. Due to the boundedness of 
D, A has purely point spectrum; in fact, one can build for this operator a complete 
orthonormal set of proper eigenfunctions {Fk)keK. with eigenvalues {ul)k£K, where 
/C = {1,2,3,...} and the labels are chosen so that 0 < cui ^ a ;2 ^ ^ ... (with 

the possibility that some of these inequalities are equalities, to deal with the case of 
degenerate eigenvalues). 

Let us discuss the smoothness properties of the eigenfunctions and derive some norm 
bounds for them. Clearly, for each k G {1,2,3,..} and r G R, we have G W 
and \\Fk\\r = = ujI- If r > | + j these facts and the imbedding 

(IRT^ imply Fk G since this holds for each j G N, we have 

Fk G ^“(n) . (B.22) 

To go on, let us note that fIB.lQp means the following: for each j G N and r > j + |, 
there exists a constant Aj^r G ( 0 , +oo) such that 

ll/llc. ^ A,-,||/||, for all/G-HD (B.23) 

This statement with f = Fk gives ||Tfc||(^j ^ ^ N, r > j + | or, 

equivalently, 

WFkWci ^ Air, for all j G N, 77 > 0 (B.24) 

(where stands for Ajj^d/ 2 +ri)- These results should be kept in mind in the sequel, 
together with the already mentioned Weyl asymptotics (13.101) 

(Ufc ~ C for k —)■ +CX) , 

where C := 20Fr(d/2 + l)A'^ Vol(D)-A'^ (g). 

order to derive Eq. (IB.191) one proceeds similary to footnote [32] using the embeddings 
^ H^{Q) and H'~(Al) ^ CASl); these follow from Theorem 3 on p. 155 of and some 
standard interpolation theory (see, e.g., [lOl (67] ). 

^^Concerning (13.1011 . we have already given references [32] [II]. Eq. (IB.24P is known as well in 
the literature, see [HTl l97l [T04] : here we have proposed an alternative derivation of this result just 
because it arose naturally from the general framework of the present appendix. 
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On the eigenfunction expansion for the Dirichlet kernel. In the case we are consid¬ 
ering, the eigenfunction expansion fl3.18p for the Dirichlet kernel takes the form 


H-cxd 


Ds{^,y) = Fk{x)Fk{y) ; 


k=i 


(B.25) 


hereafter we discuss the absolute convergence of this expansion with respect to the 
norm flB.lSp of x D), showing that 


+ 00 


I, ,2s| 

k=l I 


<+00 if3fis>d + |. 


(B.26) 


To this purpose, we hrst note that 


+ 00 

k=l 




2s I 


\Fk{-)Fk 


< 


\Ci 


+ CXD 

E 


+ 00 

E 


1 


l«l + l/3Kj xeO.ySO 


max max 


iiax_|a“Ffc(x)a^Ffc(y)| ^ 


k=l 


(B.27) 


max I IF; 




k\\jl\\^ k\\j2 


Let k G {1, 2, 3,...}, ji -f j 2 ^ 3 and t] > t). Due to the estimate (IB.24P we have 

II F II II F II ^ A A ! ,'^+-Jl+i2+2r; _ 

IIfc IIJ1 IIII j2 ^ ■'^2^ ’ 

but (recalling that ^ Wi) 


u 


ji+h _ 
k 





u 


ji+h _ 
1 ~ 


oji 


3-31-32 


U' 


1 


whence 

^jlV ^32V , d+j+2rj 
-^^^ OJu 

, J-J1-J2 « 

Wi 


IF, 


\3l I 


\\32 




(B.28) 


Inserting Eq. flB.28p into Eq. flB.27p we obtain the following, for all r/ > 0: 


H-OO 

E 


k=l 





(^ii+j2s£i cj'( 


max 


Aj'ir; 

j-31-32 



1 


2'^s—d—j—2ri 

^k 


From here to the end of the paragraph we assume 


Ifis > d + 


3_ . 

2 ’ 


(B.29) 


(B.30) 
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hereafter we show that the series on the left-hand side of Eq. flB.29p converges for 
a snitable rj > 0, a. fact yielding the thesis flB.26p . Indeed, dne to Eq. (IB.30p there 
is 7] > 0 snch that 


J 


3?s — — -|-j— 


d 


expressing in this way, and nsing the Weyl asymptotics (Id.lOh we get 
1 1 1 

tor k —)■ -l-cxD , 


OJu 


2^s-d-j-2ri ^(l+n)d (j{l+ri)dj^l+rf 


(B.31) 


(B.32) 


which implies convergence for the series on the left-hand side of Eq. flB.29p . 

On the eigenfunction expansions for the heat and cylinder kernels. In the type of 
conhgnration nnder analysis, the expansions fl3.47p and f|3.48p can be re-expressed, 
respectively, as 


+ CXD 

K{t]x,y) = Fk{x)Tl{y) , 

k=l 

(B.33) 

+ OD 

^(t;x,y) = ^e"‘^''Ffc(x)Ffc(y) . 
k=l 

(B.34) 


Using considerations similar to the ones of paragraph c), one shows the absolnte 
covergence of these expasions in the norm flB.lSp of C^(fl x fl), for all t > 0 and for 
each j G N: 


+ CX> 


H-oo 


J2e-^hFk{-)Fk{-)\y < +00 , < +oo . (B.35) 

k=l 


k=l 


Another result. In subsection 13.121 it is stated that, nnder the assnmptions (13.91) for 
A and considering a complete orthonormal set {Fk)k=i,2,3,... with the nsnal featnres 
for this case, the fnnction 


+ CXD 

r(t;x,y):= J]e-"*'|Ft(x)||Fi(y)| 


k=l 


(B.36) 


admits a nniform bonnd 

r(t ; X, y) ^ T(t) <-f-oo for all x, y G hi and t > 0 . (B.37) 

The proof of this statements starts from the ineqnality flB.24p with j = 0 and any 
r] > 0, this implies |Ffc(x)|, |Efc(y)| ^ Aq,, , whence 

H-CO 

T(t; X, y) ^ T(t) for all x, y G hi and t > 0 , T(t) := Ag,^ ^ . (B.38) 

k=l 

Finally, the Weyl estimates fl3.10l) ensnre T(t) < -(-cx) for each t > 0. 
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C Appendix. A ((i+l)-dimensional Green func¬ 
tion and its relation with the cylinder kernel 

Let n C R'^ be an open set, and let ^ be a strictly positive selfadjoint operator in 
(keeping into acconnt snitable bonndary conditions on We consider the 
cylinder kernel T(t;x, y) := (e“*^)(x,y) (x, y G hi, tG (0,+cxo)). 

The aim of the present appendix is to illustrate the fact mentioned in subsection 
13.111 namely, the possibility to relate T to a (d + 1) dimensional Green function. 
To this purpose we consider in the domain (9 := (0, Tcxo) x G 3 (t, x) and the 
Hilbert space we introduce therein the operator 

T* := —da + A. , (C-1) 

with suitable boundary conditions on dO := ({0} x G) U ((0, +cxd) x dVL) . More 
precisely, we assume Dirichlet boundary conditions on {0} x G and the previously 
given boundary conditions for A on (0, +cx)) x dVt. The operator V is selfadjoint; 
in the sequel we will prove that it is strictly positive. 

Let us introduce the Green function 

G(t,x;t',y) :=iP-i((t,x),(t',y)) ^ {5,5^\V-A, 5^) ; (G.2) 

this is characterized by the equation 

(-9tt + Mx)G(t,x;t',y) = (5(t-t')(5(x-y) , (G.3) 

and by the boundary conditions prescribed on dO. We claim that the cylinder kernel 
T is related to G by 

at/G(t,x;t',y)|j,^o = T(t;x,y) . (G.4) 

To prove this (and the previous statement on the strict positivity of P), let (Tfc)fceyc 
be a complete orthonormal system of eigenfunctions of A with corresponding eigen¬ 
values (cn^jfcgyc; clearly, the functions t G (0,-|-oo) i-G- sin (At) (A G (0, +cx))) 

are a complete orthonormal system in L^((0,-|-cxo)) and are eigenfunctions of —dn 
vanishing for t = 0. These facts ensure that the family of functions 

Y(A,fc)(t, x) := sin(At) Ffc(x) for (A, k) G (0, -Fcx)) x K. (G.5) 

is a complete orthonormal system of eigenfunctions of P, with 

VY(^x,k) = {X^ + ool)Y^x,k) . (C.6) 

We recall that we are assuming ojI ^ for some e > 0; the eigenvalues (A^ -t- col) 
also have as a lower bound, so V is strictly positive. 


73 




The Green function of Eq. flC.2D can be expressed via the equation 


^(hxG^y) = 


r dX dk 

{0,+oo)xK 


sin(At)Ffc(x) 


sin(At')Efc(y) , (C.7) 


which, evaluating explicitly the integral in A ([3), reduces to 

G(t,x;t',y)= [ dk 
JK 


p—ajfc(t+t') _ 

^F,(x)F,(y) 


'2,ojk 


(C.8) 


To go on, let us differentiate both sides of Eq. fIC.Sp with respect to t'; this gives 


at'G(t,x;t',y) = / dk 


IfC 


2 2 


Ffc(x)Ffc(y) . (C.9) 


Setting t' = 0 (and recalling that t > 0), the last equation yields 


at.G(t,x;t',y)|^,^o= [ dfce-‘^'=*Ffc(x)F,(y) . (C.IO) 

Jk 

The right-hand side of the above equality is just the representation of the cylinder 
kernel T given by Eq. 03.4811 : thus we have proved Eq. PC.411 . 

An example. Hereafter, the approach based on PC.411 is used to derive the expres¬ 
sion P3.60p for the cylinder kernel T in the case where 

Q-=K^, A:=-A. (C.ll) 

The associated {d + 1) dimensional domain and the operator V are 

O := (0, +oo) X , P := -d^ - A = -A^+i , (C.12) 


with Dirichlet boundary conditions on dO, which coincides with the hyperplane 
{0} X R'^; of course, in the above Ac^+i indicates the {d + l)-dimensional Laplacian. 
The Green function Q of — A^+i on the half-space O can be obtained by the familiar 
method of images from the Green function Qq of — A^^+i on the full space thus 


G{t, x; t', y) = Qoit, x; t', y) - Qo^t, x; -t', y) , 


(G.13) 


^®Just observe that, by symmetry arguments and the residue theorem, we have 

( 


sin (At) sin (Ah) 
cLX 


A2 


sin(At) sin(At') 
ci\ 


A2 


TT 




'j- 
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^ln((t-tr + (x-yn 


^o(t,x;t',y) : = 


r( 


d+l' 


{d — l)27r'*2^ ((t — t')^ + |x — yP) " 2 " 


if (i = 1 , 
ii d > 2. 


(C.14) 


Inserting Eq.s flC.lSp flC.14p into Eq. flC.4p we obtain for the cylinder kernel T the 
expression fl3.60p 


^(t;x,y) = 

(in all dimensions, including d = 1 ). 




/ fO I I I 

TT 2 + |x — y I 


d+l 


D Appendix. Derivation of Eq. (l3.90lj 


Let t h{t) be a complex-valued function, analytic in a neighborhood of [0, -|-oo) 
and exponentially vanishing for 3fJt —)■ -l-cxo. For any given s G C with 3fJs > 0, 
consider the integral 

I{s,h) ■= [ dtt^-^h{t) , (D.l) 

Jf} 

where is a Hankel contour (see below Eq. fl3.89p and Fig. [T]on page [32] for the 
description of this path); the complex power in the above equation is defined 
following Eq.s (13.151) (13.161) . Hor any h > 0, is homotopic to the path described 
as follows: 

U u , with (D.2) 

;= {t G C I t = n ± zh, n G [0, + 00 )} , 

;= {tG C I t = 9 G (7r/2,37r/2)} . 

Due to this remark and to the analyticity of h we can replace Sj with in Eq. 
(EH; SO, for any d > 0 we have 


I(s, h) = Ij (s, h) + (s, h) + 4 (s, h) , (D.3) 

^+00 ^37r/2 

/±(s, h):=T dv {v ± i6y-^h{v ± iS) , / 5 °(s, h):=i dO (h e*^)*h(h e*®) . 

Jo J-kI2 

We are now going to consider the limit 5 —?• O’*". Notice that in this limit (n-|-zh)^“^ —?■ 
while {v — iSy~^ —)■ = e 2 i 7 rs.ys-i. j^ioreover, h{v ± id) —)■ h{v). Due 

to these results, we easily infer 


lim itls, h) = 

5^0+ " 


^ + 00 


dv n® ^h{v 


(D.4) 


lim I7(s,h) = e®” 

5^0+ 


b + CxD 


dv v^ ^h{v 


(D.6) 
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Passing to the integral Ig{s,h), noting that \h{6 e^^)\ ^ C for small 6 and recalling 
that 3fJs > 0 by hypothesis, we obtain 

r-37r/2 




5Rs 


dd ^0 for 5 ^ 0+ 


/7r/2 


Summing up, in the limit h —)■ O’*" one obtains from Eq. flD.Sp that 


r‘+oo 


/(s, h) = - 1) 


dv ^h{v) ; 


(D.6) 


(D.7) 


noting that — 1 = sin(7rs), the above relation yields 


b + CxD 


/(s, h) = sin(7rs) 


dv ^h{v) , 


(D.8) 


which is equivalent to Eq. fl3.90p for 3fJs > 0, s ^ {1, 2, 3,...} . For more information 
concerning the Mellin transform and its contour integral representations see, e.g.. 


E Appendix. Derivation of Eq.s (13.1121 - 13.1141) 


We refer to the framework of subsection 13.181 about the slab hi = hli x R'^^; we retain 
all the assumptions and notations of the cited subsection. In particular, {dki)kieKi 
is a complete orthonormal system of eigenfunctions of Ai with related eigenvalues 
zul^. A complete orthonormal set of eigenfunctions {Fk)keic with eigenvalues {ojDkeic 
for the operator A in L‘^{Q) is given by 

o*k2-X2 


^fc(x) = 5'fci(Xi) 


i^k = + |k2 


(27r)d2/2 ’ 

for X = (xi, X 2 ) and k = {ki, k 2 ) G/Ci x 


(E.l) 


In the present case, the eigenfunction expansion (13.181) of the Dirichlet kernel at any 
two points X = (xi, X 2 ) and y = (yi, y 2 ) can be re-expressed as follows: 

gik2-(x2-y2) 


L>,(xi,X2;yi,y2) = 


f dki dk2 

/yCixR‘i2 (^fci + |k2p)^ 


dki (xi)5'fci(yi) 


dki dh 




5'fci(xi)5'fci(yi 


(27r)'^2 
giroj.^h-(x2-y2) 

(27r)^2 


(E.2) 


where, in the last passage, we performed the change of variables k 2 = Wk^h. On the 
other hand, it is known that, for any z G R'^^, 


dh 


o^hz 


/R'i2 


(27r)'^2 (|h|2 + l)* (27r)'^2/2 2"-ir(s) 


^,_ih(|z|) ifsGC, sRs>f, (E.3) 
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with Ky denoting the modihed Bessel function of the second kind of order z/ G C 
(see, e.g., [31801 [103]). Thus 


L>5(xi,X2;yi,y2) = 


^ ^ (^fcl|x2-y2| 

2s—di (^l) u/ci (yi) 




A' i(rat,|x2-y2|) . 


(E.4) 


(2i)*/2 2-ir(s) 

For the sake of brevity, for any z/ G C, we put 

012 : (0, +CX)) -)■ C , z ^ 0i2(^) := Ky{^/z) ; 


(E.5) 


due to the asymptotic behaviour of the Bessel function Ky near zero (see [80], p.252, 
Eq.10.30.2), for 3fJz/ > 0 this function has continuous extention to 2 ; = 0, given by 

0,(0) = 2"-ir(z/) . (E.6) 


To proceed, note that Eq. flE.4p can then be rephrased in terms of the function 0^, 
as follows: 


l),(xi,yi;g) = 


T>s(xi,X2;yi,y2) = A(xi,yi; |x2-y2 

21 -S r 




dki — r 2 ^ (E-7) 

(?fci(xi) dkAyi) Q) ■ 


This proves Eq. fl3.112p . also giving an explicit expression for the function Dg. 
To proceed note that, due to the well-known facts on the derivatives of the Bessel 
function Ky for any z/ G C (see [SU], p.252, Eq.10.29.4), we have 


d(3y 

dz 




i;=Vi 2\pZ 


= - ^ z‘^2'iF,_i(v^) = - ^ 0,-i(^) . 

Using the above identity it can be proved by induction that 

d”0 / iV 

A) = - j 0,-.(^) for n G {1, 2, 3,...} 

this fact, along with Eq. flE.Tp . implies 


(E.8) 


(E.9) 


^-(xi,yi;g)= ^ ^ 


dki 


dq^ 




ZXZ^^fci(xi)5^A:i(yi) '&g_'h_A'^k,(l) ■ (E.IO) 
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Now, recalling Eq. flE.6p we conclude that, for any n G {1,2,3,...} and any s G C 
with > j+n, 


dq^ 


(xi,yi;0) 


(-l)"r(i,-|-n) f dh , , —, 

(4ir)*/2 4"r(s) h, 


(E.ll) 


Due to a representation analogous to fl3.18p holding for the reduced Dirichlet kernel 
Eq. flE.lip implies Eq. fl3.114p . Finally, Eq. fl3.113p is just the case n = 0 of 
Eq. (IXmi) . 


F Appendix. Some results on boundary forces 


As in the hnal part of subsections I4.3l and l4.41 we work on a domain D with Dirichlet 
boundary conditions. 

F.l Derivation of Eq. (I4.25P for the pressure. We start from Eq. fl4.24p . 

holding for general boundary conditions. In the Dirichlet case that we are consid¬ 
ering, only the terms involving mixed derivatives (with respect to both x and y) 
of the Dirichlet kernel yield non-vanishing contributions on the boundary dQ 
thus, for any xG(9D, Eq. (I4.24p reduces to 


Pii^) = 


d^iyj Dji +1 (x, y) 


n-^(x) . (F.l) 


y=x 


We now claim that the terms proportional to ^ in Eq. fIF.ip vanish, i.e., that 


5ij dyi — d^iyj) Du±i (x, y) (x) = 0 for all x G 9D ; 


(F,2) 


y=x 


this will yield Eq. fl4.25p . Using the eigenfunction expansion fl3.18p for the Dirichlet 
kernel, we see that Eq. flF.2p holds if we are able to prove that 

s.,a‘F,aM-\d,F,a,T,-ld,F,my)nH^ h. keK:,^edn. (f.3) 

The simplest way to prove Eq. flF.3p is to derive the following, equivalent statement: 
for all fc G /C and all (sufficiently smooth) vector held © = (©*) : dfl —)■ R'^, 


'an 


da & I % a'FtdiFt - - diFidjFt - - djFAF^ ) rf = 0 . 


(F.4) 


39 


One can easily infer this statement using the eigenfunction expansion (13.181) . 
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Let us sketch a derivation of Eq. flF.4p . for given k & K, and © : dVL —> R'^. To this 
purpose we consider a smooth extension of © to a vector field © : dVL U R'^ 
and fix the attention on the integral 

^ / dx ^ / dx d0d^&) |F,|2 (F.5) 

^ do ^ do 


(note that dWj = A). We re-express both sides in the above identity integrating 
by parts with respect to all the derivatives appearing therein (0), considering them 
in the two orders proposed in the two sides; some of the boundary terms arising in 
this way vanish since Fk is zero on dVt. The difference between the two expressions 
thus obtained, which is obviously zero, is found to coincide with the left-hand side 
ofEq. dEl- 

F.2 Derivation of Eq. (I4.29|) . Let us stick to the framework of subsection I4.4l in 
which the domain id is bounded, and Dirichlet boundary conditions are prescribed; 
the operator A = —A -|- V, acting in L^(Q), has a complete orthonormal system of 
eigenfunctions F/. with eigenvlaues labelled by a countable set /C. We consider, 
for small e > 0, a deformation of the domain id of the form (I4.261I4.27|1 . controlled 
by a vector field © on R'^. The operator Ae '■= —A -|- V acting in has a 

complete orthonormal system of eigenfunctions with eigenvalues 
In subsection 14.41 we have already considered the regularized bulk energy correspond¬ 
ing to id^; this is (see Eq. fl4.28p i 




1 — u 
2 


keK 


We now consider the limit e —)■ 0, and expand everything to the first order in e. Eq. 
fl4.29p that we want to derive concerns the expansion in e of the bulk energy F“; as 
already mentioned, Eq. fl4.28p can be used to make contact with the expansion of 
the eigenvalues on which we now fix our attention. 

The variation of the eigenvalues under a deformation of the spatial domain for the 
Dirichlet Laplacian (or similar operators) has been the subject of classical investi¬ 
gations. Here we refer to the book of Rellich [85] (see Chapter II, at the end of §6), 
whose results can be expressed in this way with our notations: 

= ^k + ^^l + 0{e^) with wl := {Fk\BFk) , (F.6) 

where B is the selfadjoint operator in L^(fd) defined by 

Bf := di (id^&+ d^&) a,/) + Q Adi&^+ &^deV^ f (F.7) 

^°See the footnote on page 15^ 


79 



















(as matter of fact, [HS] gives the expression of B ioi V = 0, but the extension to a 
nonzero V is straightforward). Keeping in mind these facts, we return to Eq. fl4.28p 
for the regularized bulk energy; this implies 


= E“ + ee + 0(e2) , e 

To go on, we note that the dehnition of tu^ in Eq. 
yields 


2 ■ 

k&K 

fIF.bp . with B as in Eq. 


(F.8) 

(EH). 


wl= d:si Fh 


d.{(d‘6>+d’6‘)diFt) + {-Ad,6‘+e%v)Ft 


(F.9) 


The above result can be re-expressed in terms of surface integrals on dQ via suit¬ 
able integrations by parts (0); while making these computations, one must use the 
identity AF^ = (V—ojDFk (and its complex conjugate), and recall that F^ vanishes 
on dQ. In this way we obtain 


= / (ia(x) n^& 


6,,d^FAFk 


{diFk dj Fk+dj FkdiFk) 


(F.IO) 


Jdn 

We plug this relation into Eq. fIF.Sp . exchange the summation with the integration 
and use the expansion fl3.18p . which in this case reads T*s(x, y) = Yhk&K ^ -^fc(x)Ffc(y) 
in this way we infer 


= —(1—m) hF / da{x)n^{x) 6*(x) ■ 

Jan 

+ i4v)-0!i±i(x.y) 


(F.ll) 


y=x 


Now, comparing the above result with Eq. fl4.25p for the regularized pressure we see 
that 


e = -(1 -n) / da(x) 6Xx)p“(x) ; (F.12) 

Jan 

the first equality in flF.811 and Eq. flF.12H give the thesis fl4.29p . 


G 


Appendix. Derivation of Eq. 


([52T]) 


Consider the framework developed in subsection 15.21 for a fundamental operator A, 
such that <j{A) C [0, -|-oo) and 0 is in the continuous spectrum of A. Herefter we 
are using the deformed fundamental operator 

As:={VA + ef. (G.l) 

'^^See again the footnote [H] on page [351 
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We already observed in subsection 15.21 (see Eq. fl5.17|) i that the cylinder kernels T 
and T^, respectively associated to A and Ae, are related by 


= e ^‘T(t;x,y) ; 


(G.2) 


we also showed (see Eq. (15.181) 1 that, assuming the map t T(t;x, y) (for fixed 
X, y G ff) to admit an extension in a neighborhood of the real half-axis [0, -l-cxo) to 
a meromorphic function of t, we have 



(G.3) 


The integral in the right-hand side of the above equation is an analytic function of s 
on the whole complex plane, while the gamma function is meromorphic with simple 
poles at positive half-integer values of s. Taking into account these facts, hereafter 
we show how to evaluate the renormalized kernels 



(G.4) 


considering, separately, the cases Sq = —n/2 {n G {0,1,2,...}) and Sq = n/2 {n G 
{1,2,3...}). Putting together the results obtained for sq = —'^/2 (n G {0,1,2,...}) 
and for sq = 1/2, we will finally obtain the proof of Eq. fl5.2ip . 

We remark that, for sq = ±1/2, the above renormalized kernels coincide with the 
functions introduced in Eq. (15.Op . 

Case 1: Sq = , n G {0,1,2,...}. The right-hand side of Eq. (IG.3P is clearly 

an analytic function of s for < 1; thus, for n G {0,1,2,...} Zlg(x,y) has an 
analytic continuation at s = —n/2, hereafter indicated with that is 

simply obtained substituting this value of s in the integral representation (IG.3D . 
The resulting integrand is meromorphic so that, by the residue theorem. 


E>l|(x,y) = (-l)”r(n + l)Res(^r(^+i)e-"‘T(t;x,y);0) . (G.5) 


Of course, in the present case the prescription of taking the regular part in Eq. (15.Op 
is pleonastic, and the cited equation is reduced to 


^il(x,y) = lim T)l™(x,y) . 

2 e^0+ 2 


(G.6) 
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On the other hand, computation of the previous limit gives (Q) 

(x, y) = (-1)" r(n + 1) Res T(t; x, y); o) . (G.7) 

The above result can be reformulated in terms of the modified cylinder kernel T 
associated to A. Indeed, recall that T(t;x,y) = —x,y) (see Eq. fl3.54p i 
and note that, for any pair of functions f,g meromorphic near a point to, we have 
Res{fg'; to) = — Res(/' 5 f; to); these facts (and the standard identity zT{z) = r(l+; 2 )) 
give 

(x, y) = (-1)^+1 r(n + 2) Res (t-(’^+ 2 ) ; x, y); o) . (G.8) 

Case 2: Sq = , n G {1, 2, 3, ...} (with a special attention for the snbcase 

n = 1). A substantial difference occurs with respect to Gase 1. In fact, due to the 
gamma function appearing in Eq. flG.3p . the function i7^(x,y) described by this 
equation has a genuine singularity at s = n/2; in order to remove this singularity, 
it is essential to retain only the regular part in Eq. fl5.9p . To this purpose, for s 
in a neighborhood of n/2, we introduce the variable u := 2s — n and note that, for 
u —y 0 , 


2s—n ^—2ins 


K 


r(i-2s)r""^ = 


in—1 


(n —1)! [u 


+ (^ln(Kt) + 'Jem - iir - + 0{v) 


(G.9) 


where 'Jem — 0.577216 is the Euler-Mascheroni constant and, for m G {0,1,2,...}, 
■= \ (-^0 := 0) denotes the m-th harmonic number (@). 

'^^To prove this, one can proceed as follows. In the present case the cylinder kernel is assumed 
to be a meromorphic function of i, so there is an expansion r(i;x,y) = for 

some q G 7i (converging, at least, for small t; note that under the assumptions for the validity of 
(13.641) we have q = d); of course ^~k\ Cauchy formula for the product of 

two series yields 


q+n 


Res(^i ''*T(t;x,y);0) = ^^-^e 9 +„_fe(x,y) ed+„(x, y) = Res(^t fo+^^r(t; x, y); 0^ 


fc =0 


This proves Eq. (IC.7p . 

To obtain Eq. (IC.9p . one uses the following well known facts [80]: for m G {0,1, 2,...}, 


r(—u — m) = 


i-irn-u) 


{u + + m) = rn\ l + HmU + 0{u^) 


(■u+1)...(u+to) 

r(-u) = ---j + 0{u) , (Kt)" = l + (ln(/ct)-i7r)u + 0(u") for u^O. 
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It follows that 



1 r 


In(Kt) + 7 - fvr - Hn-i^ T(t; x, y) . 


(G.IO) 


According to Eq. (15.91) . the renormalized function y) is the limit £ O’*" of 

2 

the above expression. Under suitable hypotheses on the behaviour of T for 3fJl +cxd 
(namely, |T(t;x,y)| < U for some C,a>0), we can exchange the limit and 

the integral to obtain 


D 


(k) 

n 

2 


(x,y) 


2m Jsj (n-1)! 


In(Kl) + 7 - iTT - hfn-i) T(t; X, y ) ; 


(G.ll) 


the term ln(/i;t) prevents us from using the residue theorem, so we must hnd alter¬ 
native ways to evaluate explicitly the above integral. In the special case n = 1, we 
can proceed as follows. First we recall that T(t;x, y) = —cltE(t;x, y) and integrate 
by parts Eq. flG.lip to obtain 


D("Vx,y) = - 

i ^ 27rz 


dii T(t;x,y) ; 


(G.12) 


the resulting integrand is meromorphic in t so that we can resort to the residue 
theorem to obtain 

L>l''^(x,y) = Resfr^T(t;x,y) ;0) . (G.13) 


Conclusion. Eq. fIG.Sp for with n G {0,1,2,...} and Eq. (IG.lOp for D'l 
prove Eq. (I5.2ip for with n e {—1, 0,1, 2,...} . 


(k) 

1 

2 
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